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Newly formed black holes are expected to emit characteristic radiation in the form of quasi-normal 
modes, called ringdown waves, with discrete frequencies. LISA should be able to detect the ringdown 
waves emitted by oscillating supermassive black holes throughout the observable Universe. We 
develop a multi-mode formalism, applicable to any interferometric detectors, for detecting ringdown 
signals, for estimating black hole parameters from those signals, and for testing the no-hair theorem 
of general relativity. Focusing on LISA, we use current models of its sensitivity to compute the 
expected signal-to-noise ratio for ringdown events, the relative parameter estimation accuracy, and 
the resolvability of different modes. We also discuss the extent to which uncertainties on physical 
parameters, such as the black hole spin and the energy emitted in each mode, will affect our ability 
to do black hole spectroscopy. 
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I. INTRODUCTION 

The Laser Interferometer Space Antenna {LISA) is being designed to observe gravitational waves in the low- 
frequency regime, between 10“® and 10“^ Hz. A leading candidate source of detectable waves is the inspiral and 
merger of pairs of supermassive black holes (SMBHs). The signal should comprise three pieces: an inspiral waveform, 
a merger waveform and a ringdown waveform. The inspiral waveform, originating from that part of the decaying orbit 
leading up to the innermost stable orbit, has been analyzed using post-Newtonian theory and black-hole perturbation 
theory, and extensive studies of the detectability of this phase of the signal have been carried out (see eg. nil and 
references therein). The nature of the merger waveform is largely unknown at present, and is the subject of work in 
numerical relativity. 

The ringdown waveform originates from the distorted final black hole, and consists of a superposition of quasi¬ 
normal modes (QNMs). Each mode has a complex frequency, whose real part is the oscillation frequency and whose 
imaginary part is the inverse of the damping time, that is uniquely determined by the mass M and angular momentum 
J of the black hole. The amplitudes and phases of the various modes are determined by the specific process that 
formed the final hole. 

The uniqueness of the modes’ frequencies and damping times is directly related to the “no hair” theorem of general 
relativistic black holes, and thus a reliable detection and accurate identification of QNMs could provide the “smoking 
gun” for black holes and an important test of general relativity in the strong-held regime Q . 

In a pioneering analysis, Flanagan and Hughes (0, henceforth EH) showed that, independently of uncertainties 
in the black hole spin and in the relative efficiency of radiation into various modes, the signal-to-noise ratio (SNR) 
for black hole ringdown could be comparable to that for binary inspiral. Consequently, both ringdown and inspiral 
radiation from SMBHs should be sufficiently strong relative to the proposed sensitivity of LISA that they may both 
be detectable with high SNR throughout the observable universe. With high SNR comes high accuracy, and hence 
the potential to measure ringdown QNMs and to test general relativity. 

The EH analysis provided some insight into the issue of detectability of ringdown waves. However, to our knowledge, 
the problem of parameter estimation from black hole ringdown with LISA has not been discussed in depth in the 
literature to date. Most existiM studies have referred specifically to high-frequency ringdown sources and Earth-based 
interferometers |E IE 0 B S ■ 

The main purpose of this paper is to discuss the measurability of ringdown QNM frequencies using LISA by carefully 
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FIG. 1: Value of Erd required to detect the fundamental mode with Z = m = 2 (detection being defined by a SNR of 10) at 
Dl = 3 Gpc. For illustrative purposes here we pick the fundamental mode with I = m = 2, but the dependence on (n, I, m) is 
very weak. The three curves correspond to j = 0 (solid), j = 0.8 (dashed) and j = 0.98 (dot-dashed), where j — j'^/M = a/M 
is the dimensionless angular momentum parameter of the hole. The “pessimistic” prediction from numerical simulations of 
head-on collisions is Erd = 10“® (as marked by the dashed horizontal line), so we should be able to see all equal-mass mergers 
with a final black hole mass larger than about ~ 10® Mq (the vertical line is just a guide to the eye). The dip in the curves is 
a consequence of white-dwarf confusion noise in the LISA noise curve. 


developing a framework for analyzing QNM radiation, and then applying the standard “Fisher matrix” formalism for 
parameter estimation |^. We will treat both single-mode and multi-mode cases. 

From a mathematic al p oint of view the excitation amp litude of QNMs is an ill-defined concept, because QNMs 
are not complete 0, MM Following IITL we will associate with each QNM an “excitation 

coefficient” that quantihes in a pragmatic (as opposed to rigorous) way the response of a black hole to perturbations 
with some given angular dependence. We will also dehne useful energy coefficients to characterize the energy deposited 
into various QNMs. Unfortunately, there is only sketchy information at present from numerical and perturbative 
simulations of distorted black holes, gravitational collapse or head-on collisions of two black holes as to what might 
be expected for the amplitudes, phases or energies of QNMs. It is clear that the QNM content of the waveform will 
depend strongly on the initial conditions and on the details of the distortion. In the absence of such information we 
will consider appropriate ranges of energy coefficients, and ranges of relative QNM amplitudes and phases as a way 
to assess the measurability of ringdown modes in some generality. 

Although we will consider a range of SMBH masses from 10® to 10® M©, we note that, for masses smaller than 
about 10® Mq, the damping time of the waves may be shorter than the light-travel time along the LISA arms, and as 
a consequence the number of observable oscillations will be so small that a Fisher matrix approach may not be fully 
reliable. We will also consider the full range of SMBH angular momenta, from zero to near extremal. 

One of our conclusions is that the prospects for detection of ringdown radiation by LISA are quite encouraging. 
Figure Q shows the value of the fraction of “ringdown energy” Crd (defined as the fraction of the black hole mass 
radiated in ringdown gravitational waves) deposited in the fundamental “bar mode” with I = m = 2 (assuming that 
mode dominates) that is required for the mode to be detectable by LISA with a SNR of 10 from a distance of 3 
Gpc. Three values of the dimensionless angular momentum parameter j = J/M"^ = a/M are shown: zero, 0.8 (an 
astrophysically interesting value), and 0.98. Recall that 0 < j < 1, spanning the range from Schwarzschild to extremal 
Kerr black holes. For SMBH masses between 10® and 10^ Mq, deposition energies as small as 10“^ of the mass should 
be detectable. We show that this conclusion is not strongly dependent on or on the overtone index n. 

We also find that accurate measurements of SMBH mass and angular momentum may be possible. For detection 
of the fundamental I = m = 2 bar mode, for example, Figure |21 shows the estimated error (multiplied by the SNR, 
p) in measuring the SMBH mass M, angular momentum parameter j, QNM amplitude A, and phase (j) (see Sec. IHII 
for definitions). For example for an energy deposition of 10“'*M into the fundamental mode of a 10® Mq SMBH with 
j = 0.8 at 3 Gpc {p ~ 200), M and j could be measured to levels of a percent; if the energy deposition is only 10“®, 
they could still be measured to 10 percent. Generalizing to multi-mode detection (specifically to detection of two 
modes with a range of relative amplitudes), we find similar results. 

In order to test the no-hair theorem, it is necessary (though not sufficient) to resolve two QNMs Q; roughly 
speaking, one mode is used to measure M and j, the other to test consistency with the GR prediction. Using a simple 
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FIG. 2: Errors (multiplied by the signal-to-noise ratio p) in measurements of different parameters for the fundamental I = m = 2 
mode as functions of the angular momentum parameter j. Solid (black) lines give paj, dashed (red) lines pam/M, dot-dashed 
(green) lines pa a/ A, dot-dot-dashed (blue) lines pa^, where au denotes the estimated rms error for variable k, M denotes the 
mass of the black hole, and A and (f> denote the amplitude and phase of the wave. 



j 

FIG. 3: “Critical” SNR pcrit required to resolve the fundamental mode (n = 0) from the first overtone (n' = 1) with the same 
angular dependence {I = I', m = m'). We assume the amplitude of the overtone is one tenth that of the fundamental mode. 
Solid lines refer to m = 1,.., 1 (bottom to top), the dotted line to m = 0, and dashed lines to m = —1,.., —I (bottom to top). 


extension of the Rayleigh criterion for resolving sinusoids, we estimate the SNR required to resolve the frequencies 
and/or the damping times of various pairs of modes, as a function of the angular momentum parameter j. For 
example, to resolve the fundamental (n = 0) mode for a given (l,m) from the first overtone (n = 1) for the same 
the critical SNR required to resolve either frequency or damping time is shown in Fig. |3 while the SNR 
required to resolve both is shown in Fig. 0] These values assume that the amplitude of the overtone is 1/10 that of 
the fundamental mode. Comparing Figs. Eland El with Fig. [7|in Sec. IIII Bl we infer that tests of the no-hair theorem 
should be feasible, even under rather pessimistic assumptions about the ringdown efficiency as long as the first 
overtone radiates a fraction ~ 10“^ of the total ringdown energy. However, resolving both frequencies and damping 
times typically requires a SNR greater than about 10^. This is only possible under rather optimistic assumptions 
about the radiative efficiency, and it can be impossible if the dominant mode has I = m = 2 and the black hole is 
rapidly spinning (solid black line in the left panel of Fig. EJ- Requiring SNRs at least as large as 10^ implies that 
resolving QNMs will be impossible for redshifts larger than about 10. 

The remainder of this paper provides details. In Sec. m we introduce our notation and formalism, and clarify 
some conceptual issues related to the QNM decomposition of gravitational waveforms. In Sec. cni we compute 
the single-mode SNR assuming that only one mode dominates the ringdown. As a first step we update the FH 
analysis of LISA's SNR for detection of single-mode waveforms. With respect to FH we use a better semianalytic 
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FIG. 4: “Critical” SNR pboth required to resolve both the frequency and the damping time of the fundamental mode (n = 0) 
from the first overtone {n' — 1) with the same angular dependence (Z = I', m = m'). We assume the amplitude of the overtone 
is one tenth that of the fundamental mode. Solid lines refer to m = I ,1 (top to bottom), the dotted line to m = 0, and dashed 
lines to m = —1,.., —I (top to bottom, unless indicated). In the color versions, we used black for the modes with I = |m|, red 
for those with 0 < \m\ < I and blue for m = 0. For I = m the critical SNR grows monotonically as j —> 1. 


approximation of the LISA noise curve, we consider different QNM frequencies, and we compare with the expected 
SNR for inspiral as computed in 0]. We explore the angular momentum dependence by considering black holes with 
j = a/M — 0, 0.8, 0.98, and confirm the FH expectation that angular momentum does not have a big effect on the 
SNR. Uncertainties in the ringdown efficiency Crd have a larger impact, since p ^ CYl we assess the 

accuracy of parameter estimation in single-mode situations, revisiting the analyses in [a la 05 show that it is 
in general very good. A more detailed analysis shows that, for counterrotating (m < 0) and axisymmetric (m = 0) 
modes, rotation doesn’t necessarily help parameter estimation, and for some counterrotating modes the error can even 
blow up at “critical” values of the black hole angular momentum. Following preliminary remarks in Sec. 0 describing 
mathematical issues in and model predictions for multi-mode ringdowns, in Sec. EH we generalize to a two-mode 
analysis, computing the SNR and errors on parameter estimation. We find that it is computationally convenient to 
treat separately the case of waveforms in which the modes have different angular dependence {I ^ V oi m ^ m!) and 
the case where the overtones have the same angular dependence as the fundamental mode. In Sec. Em we determine 
the minimum signal-to-noise ratio required to discriminate between different quasi-normal mode pairs. Conclusions 
and perspectives for future work are presented in Sec. Emi 

In the Appendices we collect various numerical results and technical calculations. Appendix^discusses the accuracy 
of a useful semianalytic approach to the calculation of the SNR (the d-function approximation, valid for ringdown 
signals with large quality factor). In Appendix IBl we present an explicit calculation of the Fisher matrix using 
different conventions and different numbers of free parameters in the ringdown waveform. Appendix o describes 
our semianalytical model of the LISA noise. Appendix IDI discusses a particular aspect of our criterion for resolving 
normal modes. Finally, AppendixEllists for reference the values of the complex frequencies and the angular separation 
constants of the spin-weighted spheroidal harmonics for selected normal modes, and also displays analytical fits of the 
QNM frequencies accurate within a few percent. 


II. BLACK HOLE RINGING: PRELIMINARIES 
A. Optimal mass range for ringdown detection by LISA 

During the ringdown phase, perturbations of a Kerr black hole die away as exponentially damped sinusoids, whose 
frequencies and damping times are given by (complex) QNM frequencies |^ . We decompose the perturbations in 
spheroidal harmonics Sim{t,(3) of “spin weight” 2 |^l^ . where I and m are indices analogous to those for standard 
spherical harmonics, and i and (3 are angular variables such that the azimuthal, or f3 dependence goes like e'^. For 
each (Z, m) there is an infinity of resonant quasi-normal frequencies |2n|. which control the intermediate time behavior 
of the signal. We label each of these frequencies by an overtone index n such that the mode with n = 0 has the longest 
damping time, followed by n = 1 and so on. Thus, in the end QNM frequencies are parameterized by three numbers: 
I , m and n. Now, the time dependence of the signal during ringdown is of the form but since u = uJimn + i/Rmn 
is in general a complex number, we will follow the usual conventions and write this as sin (wjmnZ^ + 

or cos (ta/mni + where oJimn = ‘2'n'fimn is the mode’s real part and Timn is the damping time of the 
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FIG. 5: Frequency fimn (left) and qnality factor Qimn (right) for the fnndamental modes with / = 2, 3, 4 and different values 
of m. Solid lines refer to m = Z,1 (from top to bottom), the dotted line to m = 0, and dashed lines refer to m = —1 ,—I 
(from top to bottom). Quality factors for the higher overtones are lower than the ones we display here. 


oscillation. We will also define the quality factor of a QNM as 


Qlmn — '^flmn'^lmn — ^lmn'^lmn/‘^ • 


( 2 . 1 ) 


In our analysis of the detectability of ringdown radiation, we will assume that the signal lasts for at least one light- 
propagation time corresponding to the LISA arm length L ~ 5 • 10® m, or Tught = Ljc^ 16.68 s (shorter-lived signals 
may require specialized detection techniques). This places a rough lower limit on the black hole masses that are 
relevant. To see this, we note that the fundamental mode of a Schwarzschild black hole corresponds to an axially 
symmetric (m = 0), quadrupolar (I = 2) perturbation with frequency 

/200 = ±1.207 • 10-2(106Mq/M) Hz , (2.2) 


and damping time 


T 200 = 55.37(M/1O®M0) s. (2.3) 

For rotating holes, the dimensionless frequencies {MiOimn) and quality factors for the fundamental modes for 
Z = 2, 3, 4 are shown as a function of j in Fig. [3 Although the quality factors and damping times for corotating 
(m > 0) modes may be larger for rapidly rotating holes, the effects are not dramatic: for j = 0.80 the damping time 
^"220 = 65.18 (M/IO^Mq) s, and for j = 0.98 the damping time r 22 o = 127.7 (M/IO^M©) s. Accordingly we will 
restrict our attention to masses larger than 10® or a few times 10^Mq. 

We can also estimate an upper limit for masses to be considered by noting that LISA's low frequency noise 
may provide a lower cutoff at 10“'^ Hz. Equation (E3 then gives a mass upper limit of around 10®Mq; if the 
LISA performance can be extended down to 10“® Hz, then masses as large as 10^Mq may be detectable. Again, these 
rough bounds are not terribly dependent on the black hole spin or the mode. 


B. Quasinormal mode decomposition and polarization of the waveform 


The plus and cross components of the gravitational waveform emitted by a perturbed Kerr black hole can be written 
in terms of the radial Teukolsky function Rimuj as 


hj^ + ihx 




E 


Slm {l,P) Rl muj (j ") ■ 


(2.4) 


The radial Teukolsky function Rimu ~ as r ^ 00 , where is a complex amplitude [^ . 

We assume that the gravitational wave signal during the ringdown phase can be expressed as a linear superposition 
of exponentially decaying sinusoids. QNMs are known not to be a complete set, and thus such an expansion is 
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FIG. 6: Fundamental I = 2 QNM frequencies of the Kerr black hole in the range j € [0,0.99]. Solid (blue) lines correspond 
to m = 2, dot-dashed (black) lines to m = 1, dotted (red) lines to m = 0, dashed (green) lines to m = —1, and short dashed 
(violet) lines to m = —2. 


not well defined mathematically. However numerical simulations of a variety of (perturbative and non-perturbative) 
dynamical processes involving black holes show that, at intermediate times, the response of a black hole is indeed well 
described by a linear superposition of damped exponentials. We defer further discussion of the meaning of the QNM 
expansion to Sec. m For the time being, we just assume that we can write the gravitational waveform as a formal 
QNM expansion (rather than as a standard Fourier expansion in the real frequency cu) of the Teukolsky function, so 
that we can replace Eq. 12.411 by 

= J E (b /3) , (2.5) 

Imn 

where n denotes the overtone index and from here on the coordinate t stands for the retarded time t — r. We write 
the complex wave amplitude in terms of a real amplitude Aimn and a real phase (j)imn, and to follow the FH 

convention we factor out the black hole mass M: Zf^^ — M■ In this way we get 

h+ + i/lx = E ^ _ (2.6) 

Imn 

In this expansion the spheroidal functions Simn = Simiauumn) are evaluated at the (complex) QNM frequencies, so 
they are complex numbers (henceforth we drop the (i,/3) angular dependence on the Simn)- 

One frequently finds in the literature the astrophysically reasonable assumption that only the I = m = 2 mode is 
present in the waveform. This viewpoint has two conceptual flaws. 

First, QNMs of Schwarzschild and Kerr black holes always come “in pairs”. In the Kerr case, for a given {l,m) and 
a given value of a = jM the eigenvalue problem admits two solutions: one with positive (real part of) the frequency, 
the other with negative real part of the frequency and different damping time. To illustrate this property, in Fig.|^ 
we show the fundamental Kerr QNM with I = 2 and different values of m. Positive-m frequencies are related to 
negative-m frequencies by a simple symmetry property: one can easily see from the perturbation equations that to 
any QNM frequency characterized by {l,m,n) there corresponds a QNM frequency characterized by such 

that 


^Imn — ^l — mn : ^— j -^Imn — Al — jnn : (2-7) 

(here Aimn is the angular separation constant, not to be confused with the mode amplitude Aimn)- In this sense, 
any solution with positive m is nothing but the “mirror image” of a solution with negative real part and negative m. 
For m = 0 (and, of course, in the Schwarzschild case) the two “mirror solutions” are degenerate in modulus of the 
frequency and damping time. However, in general, a “mode with a given will always contain a superposition 

of two different damped exponentials. One of these exponentials could be invisible in the actual waveform because its 
damping time is shorter, or perhaps because it is less excited in the given physical situation, but formally one can 
never have anything like an isolated “I = m = 2 frequency” with a positive real part. This excitation of both modes 
is actually observed in time-evolutions of perturbative fields in Kerr backgrounds 1^ . 
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Second, a single-mode expansion automatically restricts attention to circularly polarized gravitational waves; more 
generally one cannot specify the polarization state of the waveform by assuming that it is described by a single QNM 
frequency. This problem has been overlooked in all previous treatments of the gravitational radiation emitted during 
ringdown. This omission has no serious consequences for nonrotating holes, but it is conceptually inconsistent when 
j ^ 0. Consider, for example, the starting point of the analysis in FH. They assume that the waveform can be written 
as [Eq. (3.15) in FH] 


, ( 2 . 8 ) 

r 

with I = m = 2. This is not a general assumption: it implies that the waves are circularly polarized. 

On the contrary, the polarization of the ringdown waveform depends on the physical process generating the distortion 
of the black hole. In fact, most studies dealing with point particles in the vicinities of black holes show that the 
component is extremely hard to excite |^ . The only exception is provided by particles in circular motion, which 
do not resonantly excite QNMs anyway. To excite QNMs one needs an object passing through the potential barrier 
peak, for instance during the merger phase. In this case the motion is almost radial, so one can presumably lean on 
the point particle results and assume hx ^ 0 (see however |23,l23 discussion of possible circular polarizations). 

Accordingly, a consistent approach to ringdown waveforms begins with a general superposition of modes, including 
the “twin” modes with frequency = —coi-mn and a different damping = Ti-mn- Then, using the symmetry 
property we can easily see that: 


hj^ + i/ix 


Imn ^ 


M 

r 


Imn ^ 


Imn ^ 


M 

r 


Imn ^ 


Si -I- a; pi(“i'mn*+'A!mn)p-*/Amn C' I 

^Imn \ '^Imn^ ^ ^Imn j 

Si -l-Al pi(-“i-mnt-|-0Ln)p-‘/Ti-mn C* I 

^Imn \ '^Imn^ ^ ^l — mn | 

Q, \ A' C* [ 

Jlmn W c ^Imn I 

Si A- A', C* I 

^Imn W •^Imn I ■ 


(2.9) 


In going from the second to the third line we relabeled m —>■ —m in the second term, and in going from the third to 
the fourth line we changed the labeling of the (arbitrary) constants, replacing by and by (t’lmn- 

So the general waveform depends on four arbitrary, real constants: Aimm 4>imn and (ji'imn each {I, m, n). 

Thus it is clear that only by combining positive and negative values of m can we require the waveform to have any 
given polarization state. In particular, if Simn = ^imn waveform becomes pure real (we 

have a “plus” state); if instead Aimne-^^^'^’^"Simn = ^imn becomes pure imaginary (we have a “cross” 

state). In our single-mode analysis we will usually write the (real) plus and cross components measured at the detector 
as damped sinusoids, specifying arbitrarily their amplitude and relative phase. More rigorously, when we write the 
waveform as a damped sinusoid we really mean that we have performed a sum of the appropriate QNM components, 
as described above. 


C. Including cosmological redshift 

The general waveform (HH) is written in the rest frame of the black hole, and thus all the quantities appearing 
there (M, ujimn and Timn) are measured in that frame. However, because of cosmological effects, in the detector’s 
frame all dimensionful quantities should be interpreted as redshifted. The prescription to include cosmological effects 
is very simple r should be replaced by the luminosity distance Dl{z), and all quantities with dimensions 

[mass]^ should enter the waveforms at the detector multiplied by the factor (1 -|- z)p. So, whenever the source is at 
cosmological distance, our r should be replaced by Dl{z), M by the redshifted mass {l + z)M^, fimn by the redshifted 
frequency (1 + z), and rimn by (1 -I- (where all quantities marked by a superscript 0 are measured in the 

source frame). In our numerical work, we use the values of cosmological parameters reported in |^. 



III. SIGNAL-TO-NOISE RATIO FOR A SINGLE-MODE WAVEFORM 


A. Analytic Results 


We begin by studying the SNR for detection of a single QNM. From Eq. we can express the (real) waveform 

measured at the detector as a linear superposition of /i+ and /ix, where, for the given mode {I, m, n), 


h 

h 


+ 


X 


r 


'^Imn^ 


^Imn (^; P') 

AX JiuJlm 

'^Imn^ 


^ Slmni^f P') 


(3.1a) 

(3.1b) 


where and 4>l^ are real, and are related to the quantities Aimn, (t^imn, and ^q. (12.91) by 

± , where the +(—) signs correspond to the -l-(x) polarizations, respectively. 

The waveform measured at a detector is given by 


h = h+F+{6s,4>sAs) + hxFx{0s,4>sAs), 


(3.2) 


where -F+ x are pattern functions that depend on the orientation of the detector and the direction of the source, given 
by 


F+{9s,4>sAs) = -{I + cos^ 9s) cos2<j)s cos2tps — cos9ssm2(j)ssin2tps , (3.3a) 

FxiSsAsAs) = ^(1+ cos^ 6 ( 5 ) cos2(/)s sin2 '!/'s + cos 05 sin2(/)s cos2V’s . (3.3b) 

To compute the SNR we will usually follow the prescription described in Appendix A of FH (henceforth the FH 
convention or FH doubling prescription) as follows: (1) Assume that the waveform for t < 0 is identical to the 
waveform for t > 0 except for the sign of t/rimn, i-C. that we replace the decay factor .virith cFAHimn_ ^ 2 ) 

Compute the SNR using the standard expression. 


= 4 


h*if)hif) 

Shif) 


df, 


(3.4) 


where h{f) is the Fourier transform of the waveform, and Sh{f) is the noise spectral density of the detector. (3) 
Divide by a correction factor of in amplitude to compensate for the doubling-up in step (1). 

In calculating the SNR, we will average over source directions and over detector and black-hole orientations, making 
use of the angle averages: {F^) = {F^) = 1/5, {F+Fx) = 0, and {\Simn\'^) = l/ 47 r. This simple averaging is feasible 
because the mode damping time is short compared to the orbital period of LISA. 

Sometimes, for comparison, we will not follow the three steps we just described, but will calculate the Fourier 
transform of the waveform by integrating only over the range t > 0. Since this was the method used by Echeverria 
and Einn [^, we will refer to this procedure as the Echeverria-Finn (EF) convention. 

In the rest of this Section we will follow the FH prescription. The Fourier transform of the waveform can be 
computed using the elementary relation 


f e‘“* (e±‘ 
J — 00 ' 




dt = 


2.lTir. 


(l/T-imn)^ + (pj±UJl mn)^ 
Then the Fourier transforms of the plus and cross components become: 


= 2b±. 


r M , 

h+ = —A 


Imn 




■M ,x 

= -i—A"" 


imn 




Imn'" 


(3.5) 


(3.6a) 

(3.6b) 
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We can directly plug these Fourier transforms into the definition El of the SNR to get 




(mV /■°° df 
\r)J, Snif) 

^ I (^+ + ^-) [■^hnn^+ + ~ ‘^■^hnn-^fmn^+^X ~ '/^Zmn)] \^lmn\ 


+26+6_ 


5R 


j+2 p2 _ j,x 2 ^2 o 2 


Imn^ X'• 


^Imn 


_|_ 2 d+ P, P^c> c 2 \ 


(3.7) 


The terms proportional to cannot be angle-averaged analytically in the usual way, so to deal with this general 
expression one must perform a Monte Carlo simulation. Given randomly generated values of the angles we can 
compute numerically the spin-weighted spheroidal harmonics at the QNM frequencies, plug the harmonics into the 
integrals, and finally average the resulting SNRs. We leave this for future work. 

However, especially for slowly-damped modes with m > 0, the imaginary part of Simn is typically smaller than the 
real part: 


Slmn ^ > %Slmn) ■ (3.8) 

We give a quantitative discussion of the validity of this approximation elsewhere . We will henceforth assume that 
the Slmn are real, so that Sf^^^ = \Sf^^\, and we can complete the angle averaging analytically to obtain 

^ lOir \ r ) Jq Shif) ^ 

X {{bl + bl) [At^^ + + 26+6_ [Atl cos(2<^+„ J - cos(2,^L J] } . (3.9) 

We expect that the resulting SNR should be reasonably close to the true angle-averaged result as long as the imaginary 
part of the harmonics is not too large. 

We make the further approximation that the damping time is sufficiently long that the frequency-dependent func¬ 
tions 5^ -I- 6^ and 6+6_ may be replaced by suitable (^-functions, namely in the large Qimn or large Timn limit, 

b^ +b‘^ ^ [S{f — flmn) + < 5 (/ + flmn)] , 

b^b_ ^ / m - flmn) + S{f + flmn)] , (3.10) 


where the normalizations are obtained by integrating over positive frequencies only. This approximation is mathe¬ 
matically, though not physically equivalent to assuming that the noise density Sh{f) is strictly constant. We then 
obtain the angle-averaged SNR, 


P^ = 


Qln 


407r^/;mn(l + 4 :Qfmn)^h{flmn) 


X 


-b 





[1 -b cos{2(j)+^J + 4Q;An] 

[1 - C0S(2^,^^„) -b AQjmn] I 


(3.11) 


For simplicity, FH also make an assumption about the relative amplitudes and phases of the waves, taking a pure 
cosine for the -b-polarization {(f'^mn = 0 ), a pure sine for the x-polarization {4>fmn — ^"^d assuming = Amn = 

Aimn- With these assumptions, the SNR takes the form, 


Pfh 


[mV r df f _^_ 

[ r J Jo Shif) 1 [(/ -b fimn)'^ + (27rTi™„)-2]^ 

(W\ QlrnnAfmn 

\ r J 20Tr‘^fimnShiflmn) ’ 


[(/ - flmn)'^ + {2TTTimn) 


(3.12a) 


(3.12b) 


where the second expression corresponds to the 5-function limit. 
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It is now useful, following FH, to define an energy spectrum through the relation 


P^ = 


1 


dE 


57r2r2 pShU) df 


df. 


(3.13) 


From Eq. (mill . we obtain 
dE ttM^P 

W ~ 

We then define the “radiation efficiency” Crd, by 

Erd = 


{{bl + b^_) + All] + 2A&- [ALl cos( 2 <^+„ J - All cos(2</.A J] } • (3-14) 


Egw 

M 


M 


dE 

W 


df. 


(3.15) 


Substituting Eq. ijnijl into 13.1511 and integrating, and comparing the result with Eq. (linrii . we find a relation 
between SNR and radiation efficiency for a given mode, in the J-function or constant-noise limit. 


Pfh = 


1/2 


/ CrdM 


mn ) / \/l -|- ^Ql. 


1/2 


‘^Qln 


(3.16) 


independently of any condition on the relative amplitudes or phases. 

With the FH choice of phases and amplitudes, the resulting energy spectrum is their formula (3.18): 


(dE\ 

1 \ 

r 1 1 1 1 

\df) 

FH 327r3A^„ 1 

1 [(/ + flmny + (27rT;m„)-2]^ [(/ - fimnY + (27rr;m„)-2]^ j 


mn Mpi 

71 


-b{f - flmn) ■ 


(3.17) 


where the second expression corresponds to the i5-function limit. Integrating the FH energy spectrum lITTTIl explicitly, 
we find that the amplitude is related to e^d by 


•^Imn — 


^^Qlmn^rd 


M fimnd + \MQimnfl 


8 erd 


(3.18) 


where the second expression corresponds to the d-function limit. 

Using our general spectrum we can relate the polarization-phase dependent amplitude to an efficiency per 

polarization by 


A~^’^ 

'^Imn 




+,'X 

^rd 


16e 


-e.x 

rd 






± cos(2<).l;0)] 


MQi 

mnfln 


(3.19) 


where the upper and lower signs refer to the -I- and x polarizations, respectively, and where the last step again 
corresponds to the i5-function limit. 

The expressions used in this Section are valid for any interferometric detector. In all of our LISA calculations 
we take into account the fact that the LISA arms form an angle of 60 degrees; as a result, when integrating our 
results with the LISA noise curve, we must multiply all amplitudes by a geometrical correction factor so that 

^y^^ = y3/2x^+.x. 

We now combine the expression 113.1611 with an analytic approximation for the LISA noise curve (here, for simplicity, 
we exclude white-dwarf confusion noise; see Appendix EJ of the form 


sr\f) 


LO-52 (JLS) 


-4 


-F 1.59 X 10-41 


Hz 


-1 


(3.20) 


Rescaling frequencies in terms of the dimensionless frequency Amn = Mujimn, and inserting redshift factors suitably, 
we obtain 


5.1 X 103 ^^y/2 ( {l + z)M \^/^ ( I Gvc \ ( So \ 2Qimn 
Eimn V0.03J V J \DL{z))\Sh{fi^p) yi + 4QL„ ’ 
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where 

Shiflmn) 5.4x10-5 fil + z)MY , . , f y 

5o V io^mq ; + 


(3.22) 


The dimensionless, mode-dependent quantities Timn and Qimn are of order unity and vary relatively weakly from 
mode to mode; for low order modes they can be determined from the analytic fits discussed in Appendix El 

So far we have confined attention to the FH convention. If we follow the alternative EF convention of keeping the 
waveform only for t > 0, and use the (^-function limit, we get the SNR 

Pef = ( 1 ^) = 40,y.J(‘r+4Ql^JS {(^^) [' + '“(2-*;„„)-20,™nSin(2CJ + 40U 

+ [l-cos(2</2r_) + 2Q,™„sin(2</))<_) + 4QLJ M3.23) 


where the additional phase-dependent term comes from the lack of time symmetry imposed on the waveform. The 
rest of the formulae in this section can be recast simply using this convention. 


B. Numerical results 

We first compute SNRs for both inspiral and ringdown for events at = 3 Gpc, corresponding to a redshift 
z ~ 0.54, based on our choice for cosmological parameters To compute the inspiral SNR we adopt the method 
discussed in Ref. P|. We perform an angle-average over pattern functions, assuming that we observe the last year of 
inspiral and that we can extrapolate the LISA noise curve down to a frequency / ~ lO”® Hz. Following the common 
practice, we truncate the signal-to-noise ratio integral, Eq. (lOl . using an upper cutoff frequency determined by the 
conventional Schwarzschild ISCO for a black hole of mass M. 

We compute the ringdown SNR for the fundamental mode with I = m = 2; calculations for different values of I 
and m yield similar results, the SNR depending mainly on the ringdown efficiency in a given mode. We use the FH 
SNR 13.1211 and adopt the ^-function approximation I3.12bll . Performing the “full” integral over the Breit-Wigner 
distribution 13.12all we obtain essentially indistinguishable results, except for a small (< 10%) disagreement in the 
mass/frequency region dominated by the white-dwarf confusion noise (this statement is made more quantitative in 
Fig- 1121 of Appendix IXll. 

The results are shown in Fig. 0 These plots can be viewed as an updated version of Fig. 6 in FH. Compared to 
FH we use a better model of the LISA noise curve (cf. AppendixEJ. In particular, a comparison of the left and right 
panels illustrates the effect of white-dwarf confusion noise on the expected SNR. In both panels of Fig.[71 the thick 
curve marked by “inspiral” represents the inspiral SNR for two equal-mass black holes with total mass M equal to 
the mass of the final black hole. 

The ringdown SNR in Fig. 0 is shown as sets of solid, dashed and dot-dashed lines, corresponding to the limiting 
case of a Schwarzschild black hole with j = 0, an intermediate rotation rate j = 0.8, and a near extremal rate j = 0.98, 
respectively. The intermediate value seems astrophysically quite plausible, based on the best available astrophysical 
observations and on theoretical models of merger and accretion (see, e.g. |^). FH considered only (j = 0.98). For the 
ringdown efficiency we show the optimistic value Crd = 3 % considered by FH, as well as a pessimistic value Crd = 0.1 %. 
The latter value corresponds to the present best estimates for the energy emitted in a maximally symmetric merger, 
i.e. in the head-on collision of equal-mass black holes (see and references therein). For unequal-mass mergers, FH 
suggest (interpolating between numerical and perturbative results) an energy scaling of the form (4^/M)^, where /j, 
is the reduced mass. 

The general features of the SNR curves are easy to understand. The SNR is basically proportional to the inverse 
of the noise power spectral density Sh{f)- It has a maximum in the mass range M ~ IO^Mq corresponding to the 
frequency / ~ 10“^ Hz at which LISA is most sensitive. If we include white-dwarf confusion noise (left panel) we 
observe the appearance of a dip in the SNR at masses M of the order of a few times 10 ®Mq. The black hole at 
our Galactic Center has an estimated mass M ~ 3.7 ± 0.2 x 10® Mq (see eg. H^), so an accurate modelling of the 
white-dwarf confusion noise might be very important for detection of black hole ringdown from galactic centers. In 
this paper, unless otherwise stated, we will include white-dwarf confusion noise in all of our numerical calculations. 

Fig .|7|illustrates that, even under pessimistic assumptions, the ringdown SNR is generally comparable to the inspiral 
SNR. Reducing the rotation rate does not have a dramatic effect, degrading the SNR of corotating modes by factors 
of order unity. The crucial element for detectability is the fraction of mass-energy e^d going into each mode. Note 
that Eq. 13.Hill implies that p ~ 
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FIG. 7: Comparison of the SNR for inspiral and ringdown waveforms. In the left panel we use the Barack-Cutler noise-curve; 
in the right panel we use the same noise curve, but we do not include white-dwarf confusion noise. The thick (black) line 
marked by “inspiral” is the (angle-averaged) SNR for the inspiral of two equal-mass black holes a.t =?> Gpc. The other sets 
of lines (red and blue in color versions) show the SNR for the I — m = 2 mode using the d-function approximation, assuming 
a ringdown efficiency = 3% and 0.1% respectively. Solid, dashed and dot-dashed lines correspond to j = 0 (Schwarzschild), 
j = 0.8 and j = 0.98 respectively. 



FIG. 8: Dependence of the SNR on redshift, for both inspiral (continuous lines) and ringdown (dashed and dot-dashed lines). 
We choose a ringdown efficiency trd = 1% and consider the cases j = Q, j = 0.8 and j = 0.98. From top to bottom the lines in 
each panel correspond to 2 = 0.54 (black in color versions; Dl = 3 Gpc), 2 = 1 (red), 2 = 2 (green), 2 = 5 (blue) and 2 = 10 
(purple). The dashed lines are obtained from the full integral, the dot-dashed lines use the 5-function approximation. 


FH and Ref. pointed out that, depending on the ringdown efficiency, there might be a “critical mass” at which 
black hole ringdown becomes dominant over inspiral. Assuming an efficiency Crd = 1% and a final black hole angular 
momentum j = 0.997, Ref. found that the SNR is greater in the ringdown signal for M > IO^Mq when 2 = 1. 
Their result is consistent with our Fig.[7| in addition, we find that this “critical mass” for the transition from inspiral 
to ringdown dominance depends only weakly on j, being more sensitive to the efficiency Crd- 

Ringdown efficiency plays a very important role in the ringdown SNR. Unfortunately, numerical relativity simula¬ 
tions do not provide us with reliable estimates of Crd 113 • Ref. |3(i| provides alternative, semianalytical estimates of 
the energy radiated in the plunge and ringdown phases for Earth-based detectors; unfortunately, an extrapolation of 
those results to binaries of relevance for LISA is not available. Given our ignorance of Crd, it makes sense to ask the 
following question: how much energy must be channelled into a given mode in order for it to be detectable? 

The answer to this question is provided by Fig. ^ in SecQ! where we assume p = 10 as a criterion for detectability 
and we show the fundamental mode of a Kerr black hole with I = m = 2. The result is encouraging: even in the 
pessimistic situation of a head-on collision Crd = 10“^, we can reasonably expect LISA to detect all mergers yielding 
a black hole of mass M > 10®M©. Even accounting for the fact that erd may be lower for unequal mass mergers, the 
prospects for detection are still encouraging. 

The cosmological redshift affects the detectability window and the SNR both for inspiral and for ringdown waves. 
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both through the decreasing signal strength with distance and through shifting the relevant frequencies to different 
parts of the LISA noise spectrum. Fig. 0 shows the results of numerical calculations of the SNR. For the ringdown 
signal, we pick a ’’best guess” for the efficiency Crd of 1%, intermediate between the 3% of FH and the 0.1% from 
head-on collisions. Each plot gives the SNR as a function of (mass in the source frame) for a different value of j 
(left to right, j = 0, j = 0.8 and j = 0.98). From top to bottom, the curves show sources at redshifts z = 0.54, 1, 2, 5 
and 10. The continuous lines are the inspiral SNR, the dashed lines are obtained from the full integral, the dot-dashed 
lines use the (5-function approximation. The inspiral SNR is (somewhat arbitrarily) truncated at the (large) value of 
the mass for which the starting frequency (which we pick to be one year before the ISCO, as in Q) becomes lower 
than 10“® Hz. 


IV. PARAMETER ESTIMATION BY DETECTION OF A SINGLE MODE 

A. Analytic results 

In this Section we will go beyond the issue of detectability and try to answer a different question: given the detection 
of a single QNM, what can we learn about the black hole parameters? To estimate the black hole parameters from 
ringdown waveforms, we use the standard technique of parameter estimation in matched filtering. By maximizing the 
correlation between a template gravitational waveform that depends on a set of parameters 0“ (for example, the black 
hole mass and angular momentum) and a measured signal, matched filtering provides a natural way to estimate the 
parameters of the signal and their errors. With a given noise spectral density for the detector, Sh{f), one defines the 
inner product between two signals hi (t) and /i 2 (t) by 


(/ll|/l2) = 2 



hlh2 + h*^hi 

Shif) ^ ’ 


(4.1) 


where hi(f) and /i 2 (/) are the Fourier transforms of the respective gravitational waveforms h(t). The components of 
the “Fisher matrix” Faj, are then given by 


Tab 


f dh dh\ 

[W 'WJ’ 


(4.2) 


In the limit of large SNR, if the noise is stationary and Gaussian, the probability that the gravitational-wave signal 
s{t) is characterized by a given set of values of the source parameters 0“ is 


p{e\s) 


p^°\e)exp 


1,69^^69^ 


(4.3) 


where 69°’ = 9° — 9°, and p^°^(0) represents the distribution of prior information. An estimate of the rms error, 
A9° = (((50“)^))^/^, in measuring the parameter 9° can then be calculated, in the limit of large SNR, by taking the 
square root of the diagonal elements of the inverse of the Fisher matrix. 


A9° = V^, E=r-% (4.4) 

The correlation coefficients between two parameters 9° and 9^ are given by 

Cab = . (4.5) 

We consider a waveform given by Eq. 13.111 . with the Simn assumed to be real, 

h+ = cos [27r/i™„t -f S'i™„(t, /?), (4.6a) 

hx = sin [27r/;„„t-h ()>()„„] (6,/3). (q jb) 

We also define 


^ /i~t~ = A~^ ^ A X = 4^ = 


(4.7) 
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where Nx is some numerical factor, and 


/ X 

^Imn 


A+ 

^Imn 


4^lmn ■ 


(4.8) 


Assuming that we know Nx and this waveform is dependent on four parameters (A+, (phnn^ j)i otherwise 

it depends on six parameters (A+, ^hnn^ j)- ^ popular choice for Nx is to assume that the 

distribution of the strain in the two polarizations mimics that of the inspiral phase: Nx = — 2(L • n)/[l + (L • n)^], 
where L is the orientation of the binary’s angular momentum and h is a unit vector describing the binary’s position 
in the sky. Fortunately, we will see that the errors have a very weak dependence on the number of parameters and 
on the (uncertain) value of the parameters Nx and ■ 

Assuming constant noise over the bandwidth of the signal, or taking the 5-function approximation, and using the 
FH doubling convention, we get the SNR fTTITl . In this approximation, errors and correlation coefficients can be 
computed analytically using Mathematica or Maple. The full expressions are lengthy and unenlightening, and we 
have implemented them numerically in a Fortran code. 

We first calculate the Fisher matrix in the parameter basis of (A+, fimn, Qimn), where it takes on a simpler 

form: 


where 


rA+A+ — 




1 


- A+fi„ 


^ A+Qi^ 


Nmr^Qlr, 

^QlmnQlr, 


(A+)2 

7 


(l -I- - 0) , 


A+ 


7 


2A+fimn 

7 


(l -I- ^Q'frnn “ P) > 
1 


2 AAQimn 1 + ^Q'imn 
7 (l + + 0) > 

7 


[(1 -I- ~ (1 ~ 4QLm)/3] 


9 f ’ 

Imn 

7 / 1 - 4Q 


2 

Imn 


2Qlmn V 4 + ‘^Qlmn 


7 


2/( 


2 

mn 


[(1 -I- - 0\ ) 


7 


1 


Vi mn Qi mn 4 + ‘^Qlmn 

7 4 


[(1 -I- iQ^rnnY ~ (4 — 


‘^Qlmn (4 + 4QLn) 


[(4 + 40LJ'-(1-12QLJ/3] , 


a = sin7sin20f^„-cos7sin2^7„, 
f} = siT?0cos2(j)^^^-cos^0cos2(j)+^^, 


407 r 2 /imn( 4 + 4 Q 2 ^„) ’ 


(4.9a) 

(4.9b) 

(4.9c) 

(4.9d) 

(4.9e) 

(4.9f) 

(4.9g) 

(4.9h) 

(4.9i) 

(4.9j) 

(4.10a) 

(4.10b) 

(4.10c) 


with cos'0 = 1/y^l + Nx, sin^ = Nx/+ Nx and A'^ = (A+)2(l -|- iV^) = (A+)^ -|- (A^)^. Note that, in this 

notation, pIjj = 7(1 +4Qf„„-/3). ^ ^ ... . nra 

We note that the Fisher matrix written in terms of the frequency and damping time is usually simpler |6|,l3 than 
that in terms of mass and angular momentum; however we prefer to deal directly with measurements of j and M. In 
Sec. EH to estimate QNM resolvability, we will work in terms of frequencies and quality factors. 

The transformation from the (A+, fimn, Qimn) basis to the (A+, M, j) basis is straightforward, namely, 

for any index k, 


^kM — {flmn/^)^k flrnn 7 

^kj = fLn^kfi^n + Q'lmn^kQlmr^ ^ (4-41) 

where = dfimn/dj and = dQimn/dj. 
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Converting to this basis and inverting the Fisher matrix, we find, to leading order in the errors 


O'M = 


aA+ = 


'^4’t 

^Imn 


and the correlation coefficients 


TjM 






A+4>1 


1 

Pfh 
1 

Pfh 
^/2A+ 
Pfh 
1 

Pfh 


^Qln 


1 + 


1 + 4/3 


Q'lrr 
MQimnf'imn 


fin 


1 + 


^^Imn 

ZI3 


^^Qlmn 
1 + 


1 + 4/3 


16Q2 


Imn 


1 - 


^Qlmn 


'^Qlmn 


= sgn(//,„„) X 1 - 


/■2 Q‘ 

J Imn^lmn 


^^QlmJLn 


= 


fMA+ — — 


1 

71 

1 

71 

a 


1 - 


1 - 


2Q 


Imn 

a 


2Qi 


1 - 6/3 
16(3Ln 
1 - 6/3 
16(3Ln 

7-8/3 

32Qf 

7-8/3 

32QLn 


3a 


a(10- 11/3) 


0{1/Q^), 


^V2Ql 


32V2Q'^ 


0{1IQ^), 
Oil/Q^), 

Oil/Q^). 


Imn 


(4.12a) 

(4.12b) 

(4.12c) 

(4.12d) 

(4.13a) 

(4.13b) 

(4.13c) 

(4.13d) 

(4.13e) 

(4.131) 


In calculating derivatives of the waveforms (lOll with respect to M and j (or with respect to fimn and Qimn), we 
have ignored derivatives of the spheroidal harmonics themselves. The Simn are functions of awimn = j^imn which is 
a function of j only. However, the Simn may be expanded in powers of jTimn, in the form 


Slnin — + {,j^l mn ) mYvra + 0(j^/m„)" , (4.14) 

where Y/m denotes a spin-weighted spherical harmonic, and c/um are related to Clebsch-Gordan coefficients. As a 
result, derivatives of the Simn with respect to j will be linear in derivatives of jiFimn, and because of the orthogonality 
of the spin-weighted spherical harmonics, inner products of Simn with S[^^ and of with itself will be at least 
quadratic in jTimn and its derivatives. At least for small jiFimn, we may expect these contributions to be small relative 
to the main contribution obtained by ignoring these derivatives. Nevertheless, the effect of this approximation should 

be explored further. _ 

The diagonal elements of the correlation matrix m = 1 for all i. The sgn{f[^^) in rjM comes from a \/fi^lfimn — 
\fimn\/flmn- implies that j and M have a positive correlation for corotating and axisymmetric modes (m > 0), but 
they are anticorrelated for counterrotating modes (m < 0): this is basically determined by the different sign of 
for the two classes of modes (see eg. Fig. El). 

The calculation using the EF convention proceeds along similar lines; the results are very similar, and they are 
reported in Appendix 0 The calculation for a six-dimensional Fisher matrix is straightforward, and is also relegated 
to Appendix lEI 

The large-Qimn expansions are typically accurate as long as Q'lmnIQi-mn is not very large (see Fig. [Till below, where 
this statement is made more quantitative). An analytic parametrization for a can be obtained by combining the SNR 
formula and with the QNM fits, whose coefficients are provided in Tables rVlIIl IIXI and IXI 

An important check is that the errors on M and j (and the correlation between these parameters), as predicted by 
Finn [^, agree with ours. This is not a trivial check, since Finn uses a different parametrization of the waveform. In 
particular, it is easy to check that Finn’s expressions for the errors, his Eqs. (4.20a) and (4.20b), agree with ours to 
leading order in so does Finn’s expression for the correlation coefficient, his Eq. (4.20e). The fact that only 

high-order corrections in Qfmn depend on the parametrization is a reassuring feature of the Fisher matrix calculation 
(see Appendix rBl). 
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FIG. 9: Qim.nlQ'lmn 1^6 fundamental mode with I = 2 (left) and I = 4. Solid lines refer to m > 0, dotted lines to m = 0, 
dashed lines to m < 0, and different shades (colors) denote different values of m, as indicated. Notice that this factor increases 
with j when m > 0 and decreases with j when m = —2. The factor blows up as j —> 0 for m = 0, and as j JctH for certain 
values of m < 0. This explains most qualitative features of the error plots. 


Some general comments on Eqs. and are in order. First, by combining Eqs. 14.12all and (I4.12bll with 

Eq. (IXnil . we see that the accuracy in measuring M and j can be very high under the right circumstances, namely. 


O'M 

W 
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6.8 X 10-3 X gi 

mn mi mn ^ 
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- 4 \ 1/2 
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>- 4 \ 1/2 


(.flmn) 

5h 


( Shiflmn) 


Crd 


V -s-o 


1/2 


1/2 


1 Gpc 

Dl{z) 
1 Gpc 


/10® mq m 
W + z)m) 

/10®Mo 

W + z)m) 


(4.15a) 

(4.15b) 


Notice that the measurement error is small (less than a percent), even under the very pessimistic assumption that 
a SMBH with M ^ 10® Mq radiates only a modest fraction Eqw ^ lO”'^M of its mass. The functions gimn{j) = 
Qimn/Q'imn himn{j) = {QimnfImn)/(Q'lmnfimn) depend On the particular mode we consider and on the black 
hole’s angular momentum; they are typically of order unity. For example, using the fitting relations in Appendix IeI 
we find that, for the fundamental mode with / = m = 2 of a Schwarzschild black hole, these factors take the values 
9imn{0) = 2.992, himniO) = 1.214. So, for a non-rotating black hole the error in angular momentum is slightly larger 
than the error on the mass. For a near-extremal black hole we have gimn(0.98) = 0.043, /iimn(0.98) = 0.234, and the 
error in angular momentum is now smaller than the error in the mass. We will see that this reversal of the magnitude 
of the errors for near-extremal black holes is typical of corotating modes, but does not hold true for counterrotating 
modes. 

To leading order in a large-Q/mn expansion the errors on angular momentum and mass, Eqs. I4.12all and 14.12bll . 
are proportional to Qimn/Q'imn- Fig.|2|we plot this quantity as a function of j for different modes. From the plot we 
can anticipate a few salient features. First of all, errors should decrease with rotation for corotating modes (m > 0). 
This was already pointed out in Refs. However, errors should increase with rotation for counterrotating modes 

(m < 0); even worse, at those “critical values” of j for which = 0 the errors for counterrotating modes blow up. 

Fig. 13 shows that, typically, this phenomenon is present for counterrotating modes with \m\ < 1/2. Finally, we can 
anticipate that Qimn/Q'imn (hence the error) will blow up as j ^ 0 for modes with m = 0. 


B. Numerical results 

Our expectations are validated by an explicit numerical calculation of the errors. We carry out this calculation in 
different ways, and in Fig. ^|we show that results obtained from these different methods are usually in very good 
agreement. The most reliable calculation is “fully numerical” (solid lines in Fig. IIOII . in the sense that it involves no 
semianalytical approximations. In this calculation we use the “complete” expressions for the errors obtained using 
Mathematica. For increased accuracy, for any given value of j we interpolate our numerical tables of the QNMs 
by fifth-order polynomials and evaluate “numerically” the derivatives and Q'l^^ by taking derivatives of these 
interpolating polynomials at the given j. Dashed lines investigate the accuracy of leading order Taylor expansions for 
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FIG. 10: Comparison between different approaches to compute the error (multiplied by the SNR) for the fundamental mode 
for {l,m) = (2,2) and (2,-1). Solid lines use the numerical implementation of the full expressions obtained by Mathematica, 
and a numerical calculation of the derivatives of Qimn and fimn- Dashed and dotted lines use expressions for the Fisher 
matrix to leading order in llQim,n', derivatives of Qimn and fimn are evaluated numerically for the former, and using the 
fitting functions given in Appendix^ for the latter. Black (red) lines refer to calculations of the error in j (M) using the 
FH “doubling prescription”. The plot shows that all these different calculations are in excellent agreement with each other. 
However, sometimes using the fits can produce only order-of-magnitude estimates of the errors: this happens for counterrotating 
modes, where an accurate calculation of the derivatives of Qimn and fimn is more important (dotted lines in the right panel 
deviate significantly from the “true” answer obtained by taking numerical derivatives of the QNM tables). Finally (and only in 
the right panel) we plot, in blue, results obtained using Eqs. (4.20a), (4.20b) in Finn’s paper a numerical implementation of 
the full expressions obtained by Mathematica, and a numerical calculation of the derivatives of Qimn and fimn- Finn’s formula 
would lie on top of the other lines in the left plot (for / = m = 2), but it gives a slightly different prediction for the errors on 
the counterrotating mode with I = 2, m = —1. 


large Qimn of the errors, such as Eqs. 14.1211 . and use these “local” interpolating polynomials to compute and 
Q'lmn- Finally, dotted lines use Taylor expansions of the errors and evaluate the derivatives and Q'lmn using the 
(somehow less accurate) “global” fits of the QNM tables we provide in Appendix IeI 

Overall, different choices for the doubling convention and/or for the number of parameters yield consistent results 
for the errors. In isolated, unfortunate cases the “global” fits of Appendix^ being valid in the whole range j G [0,1] 
but not very accurate for certain values of /, provide only an order-of-magnitude estimate of the errors. This happens, 
for example, when we consider counterrotating modes with I = 2 and |m| < 1/2 (see eg. the right panel of Fig. IIOII . 
It also happens for modes with m = 0, in the limit j ^ 0. The reason is that in these cases Qimn has a minimum, 
and Qimn/ Q'lmn blows up. This behavior is only captured by an accurate local fit of the QNM data (such as the 
polynomial fit we use to compute derivatives). The bottom line is that the global fits of Appendix I eI are generally 
inaccurate to compute measurement errors whenever the numerical QNM tables are such that Q'lmn — ^ for some 
value of j. 

Let us now turn to correlation coefficients. To leading order, the correlation coefficient between mass and angular 
momentum rjM = 1 for all modes. This high correlation between mass and angular momentum was first pointed 
out by Echeverria for the fundamental mode with I = m = 2. Echeverria suggested that, if we have some 
independent and more precise measurement of either the mass or the angular momentum (but not both) we could 
exploit this strong correlation to obtain an almost equally better estimate of the other parameter. This means that 
mass measurements of SMBHs (as inferred from the Keplerian orbits of the surrounding stars, for example) could 
be used in conjuction with gravitational-wave observations to provide accurate determinations of the hole’s angular 
momentum. Mass and angular momentum are also strongly correlated with the wave’s amplitude: to leading order, 
\rjA\ = ItmaI = l/v^ ~ 0.707. On the contrary, the polarization phase (/imn is very weakly correlated with the other 
parameters: the leading-order term is proportional to Q/mn when we use the FH convention, and to Q/mn when we 
use the EE convention, Eq. (ESJ. Independently of our convention and of the parametrization of the waveform, this 
small correlation implies that we can expect the phase to be irrelevant in measuring the black hole mass and angular 
momentum. 

In Fig. CH we present results from a numerical calculation of the correlation coefficient rjM (which is, again, 
independent of different choices on the doubling convention and/or on the number of parameters). 
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FIG. 11: Correlation between the mass and angular momentum parameters for the fundamental mode and different values of 
(l,m). For counterrotating modes j and M are actually anticorrelated {vjM < 0), and we plot the modulus of the correlation. 
The solid line uses the numerical implementation of the full expressions obtained by Mathematica, and a numerical calculation 
of the derivatives of Qimn and fimn- The dashed line uses an expansion to leading order in 1/Qimn of the correlation matrix, 
and a numerical calculation of the derivatives of Qi-mn and fimn- The dotted line uses again a Taylor expansion, but this time 
we evaluate the derivatives of Qimn and fimn using the fitting functions given in Appendix]^ The minimum in the correlation 
between mass and angular momentum (for I = m = 2) disappears when we use the numerical derivatives. The behavior of 
VjM, as computed numerically, is qualitatively consistent with the entry marked by fMa in Echeverria’s Table II 0, so the fake 
minimum really seems to be due to the inaccuracy of the fit. 


Fitf. 1121 shows the errors in different parameters rescaled by the SNR for different QNMs. All errors scale with the 
inverse of the SNR, ~ and all information about the detector is contained in the SNR. Therefore we plot the 
quantities {paj, pgm/M, pgaIA, which are, in some sense,“universal”: they do not depend on the specifics 

of the LISA noise curve, but only on intrinsic features of the gravitational waveform emitted by the perturbed black 
hole. 

In the plot, we use the FH convention and consider a four-dimensional correlation matrix, but results would not 
have changed much had we used the EF convention and/or a six dimensional correlation matrix. We use the numerical 
implementation of the full expressions obtained by Mathematica, and a numerical calculation of the derivatives of 
Qimn and fimn- Once again, results do not change appreciably if we use a Taylor expansion to leading order in 
^ I Qimn of the correlation matrix, and a numerical calculation of the derivatives of Qimn and fimn- Even if we use the 
QNM fits of Appendix lEl we get very similar results: only for counterrotating modes do those fits fail to reproduce 
the location of the peak we can see for {1 = 2, m = —1) and (^ = 4, m = —2). 

The general features emerging from Fig. 1121 agree with the expectations drawn from Fig. Errors on mass 
and angular momentum decrease with rotation for corotating modes (to > 0), but they increase with rotation for 
counterrotating modes (to < 0), blowing up at those “critical values” of j for which Q'lmn — (which occurs for 
counterrotating modes with \m\ < 1/2). The error on j goes to infinity as / ^ 0 when we consider modes with to = 0, 
but the error on M stays finite in this same limit. Errors on the amplitude A and phase (j)imn usually have a very 
weak dependence on j. 

Of course Fig.^ltells only part of the story, because it does not involve any information about the actual specifics of 
the LISA noise. Fig.imi shows the actual errors computed using the Barack-Cutler noise curve including white-dwarf 
confusion noise described in Appendix |21 We use the (^-function approximation to compute p, and consider different 
angular indices {l,m). For concreteness we assume that our source is located at = 3 Gpc and has a ringdown 
efficiency Crd = 3% for each mode. However, even in the worst-case equal-mass merger scenario (crd = 0.1%) errors 
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FIG. 12: Errors in measurements of different parameters for the fundamental mode with different values of {l,m), as functions 
of the angular momentum parameter j. Solid (black) lines give paj, dashed (red) lines P(jm/M, dot-dashed (green) lines P(ja/A, 
dot-dot-dashed (blue) lines po-<i>i„„- 


would only increase by a modest factor ~ 5.5. 

The angular momentum dependence can most easily be understood looking at the “universal” Fig. El As an 
example, focus on the mode with I = m = 2. From Fig. El we see that the error on j is larger than the error on M 
for j < 0.8, comparable to it for j ~ 0.8, and smaller than the error on M for j > 0.8: this is precisely what we see 
in the top-left panel of Fig. El Including information on the noise. Fig. El gives a good quantitative idea of the kind 
of accuracy we can achieve if we try to measure mass and angular momentum of black holes with LISA. Numerical 
results are in good agreement with our expectations based on Eqs. (I4.15all and (I4.15bll . Errors become unacceptably 
large only for black hole masses M < 10^Mq, but in general we can expect excellent accuracies: the measurement of 
a single ringing event can provide the mass and angular momentum of black holes with M > 5 x 10 ®Mq with errors 
smaller than one part in 10^, and (in the most optimistic cases, eg. for black holes of mass ^ 5 x 10 ®Mq) smaller 
than one part in 10®. 

To generalize our results to other values of Crd it is enough to recall that the SNR scales with the square root of 
erd, P ~ e^d^, so the errors scale like ■ Numerical simulations suggest that for the first overtone Crd should be 
smaller by a factor 10^ — 10^ when compared to the fundamental mode. This means that the error will be ~ 10 — 10^ 
times larger than the value plotted in Eig. El Accurate tests of the no-hair theorem involve the measurement of two 
QNM frequencies, so they could be possible only in a very limited mass range. This conclusion could be made even 
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FIG. 13: Errors in angular momentum aj (black) and mass ata/M (red) for: I = m = 2 (top left), Z = 2, m = 0 (top center), 
I = 2, m = —2 (top right); I — A, m = 2 (bottom left), / = 4, m = 0 (bottom center), I — A, m = —2 (bottom right), for a 
source at Dl = 3Gpc with Crd = 3%. Solid lines refer to j = 0, dashed lines to j = 0.8, dot-dashed lines to j = 0.98. Errors 
scale with p~^, hence they are proportional to 


more pessimistic when we consider interference effects in the multi-mode situation. 

All of our numerical results on the errors have been obtained assuming = 1, ^ ^tmn — 0 (values have 

been chosen to agree with FH), but the results we have presented so far can be considered more general and robust. 
In fact, the errors are almost completely independent of Nx and This is not only because these parameters 

only enter in subleading corrections, but also because the variability of the correction coefficients with Nx and 
is extremely weak. To see this, notice that the leading-order corrections to the errors on mass and angular momentum 
are given in the FH convention by (1 -|- 4/3)/16(5^^„, and in the EF convention, by (1 -I- -I- 2/3^)/8Q^mn I®®® ^dS- 

(IB2all and (IB2bll ] and use the fact that |q;| < 1 and |/3| < 1. 

In conclusion, not only are our results on the errors independent of the way we fold the waveform (EF/FH conven¬ 
tion) and independent of the way we parametrize the waveform (see Fig. IIOII : they are also largely independent of 
Nx and 


C. Bounding the black hole’s mass and angular momentum through detection of a single mode 

What kind of information can we extract from the detection of the frequency and damping time of a single QNM? 
Although we have parametrized our Fisher matrix formalism in terms of M and /, what we really measure are the 
frequency / and damping time r or quality factor Q of the ringdown wave. Unfortunately this is not sufficient to 
tell us the values of {l,m,n) corresponding to the mode detected, so we cannot determine the mass and spin of the 
black hole uniquely. The problem is that there are several values of the parameters {M, j,l,m,n) that yield the 
same frequencies and damping times. However, if we make the plausible assumption that the only modes likely to be 
detected are the first two overtones with 1 = 2, 3, 4, we can narrow the possibilities. 

Suppose that we measure a specific value of Qimn- In Tabled we list the different values of {j,l,m,n) yielding the 
same Qimn- The number of modes in our constrained set that correspond to that value ranges from one to about a 
dozen. Each of these modes then corresponds to a unique value of Timn = ‘2TrMfimn- From the measured fimm we 
can obtain a discrete list of provisional, accurately measured masses M. This list cannot be narrowed further without 
additional information, such as an estimate or bound obtained from the inspiral waveform, or the detection of an 
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Qimn 

{j, 1, m, n) 

20 

(0.988, 4, 3, 0) (0.987, 3, 3, 0) (0.976, 4,4, 0) 

19 

(0.986, 3, 3, 0) (0.986, 4, 3, 0) (0.973, 4,4, 0) 

18 

(0.984, 3, 3, 0) (0.984, 4, 3, 0) (0.970, 4,4, 0) 

17 

(0.982, 3, 3, 0) (0.982, 4, 3, 0), (0.965, 4, 4,0) 

16 

(0.979, 3, 3, 0) (0.979, 4, 3, 0) (0.960, 4,4, 0) 

15 

(0.976, 3, 3, 0) (0.975, 4, 3, 0) (0.954, 4,4, 0) 

13 

(0.989, 2, 2,0) (0.972, 4, 4, 0) (0.970, 4, 3,0) (0.946, 4, 4, 0) 

12 

(0.986, 3, 2, 0) (0.984, 2, 2, 0) (0.961, 3, 3, 0) (0.954, 4, 3, 0) (0.924, 4, 4, 0) 

11 

(0.985, 4, 2, 0) (0.981, 2, 2,0) (0.981, 3, 2, 0) (0.952, 3, 3,0) (0.943, 4, 3, 0) (0.907, 4, 4, 0) 

10 

(0.990, 4, 4,1) (0.977, 2, 2,0) (0.975, 3, 2, 0) (0.940, 3, 3,0) (0.926, 4, 3, 0) (0.884, 4, 4, 0) 

9 

(0.988, 4,4,1) (0.971, 2, 2, 0) (0.966, 3, 2, 0) (0.961, 4, 2, 0) (0.924, 3, 3, 0) (0.900, 4, 3, 0) (0.851, 4, 4, 0) 

8 

(0.984, 4,4,1) (0.962, 2, 2, 0) (0.951, 3, 2, 0) (0.932, 4, 2, 0) (0.900, 3, 3, 0) (0.861, 4, 3, 0) (0.802, 4, 4, 0) 

7 

(0.993, 2,1, 0) (0.989, 4, 3,1) (0.988, 3, 3,1) (0.978, 4, 4,1) (0.976, 4,1,0) (0.949, 2, 2, 0) (0.926, 3, 2,0) (0.875, 4, 2, 0) 

(0.863, 3, 3, 0) (0.794, 4, 3, 0) (0.724, 4,4, 0) 

6 

(0.987, 2,1, 0) (0.984, 4, 3,1) (0.984, 3, 3,1) (0.969, 4, 4,1) (0.929, 2, 2, 0) (0.897, 4,1, 0) (0.881, 3, 2, 0) 

(0.802, 3, 3,0) (0.757, 4, 2, 0) (0.671, 4, 3,0) (0.592, 4, 4, 0) 

5 

(0.990, 2, 2,1) (0.976, 3, 3,1) (0.974, 4, 3,1) (0.972, 2,1, 0) (0.954, 4,4,1) (0.892, 2, 2, 0) (0.859, 4, 0, 0) 
(0.785,3,2,0) (0.688,3,3,0) (0.620,4,1,0) (0.492,4,2,0) (0.410,4,3,0) (0.338,4,4,0) 

4 

(0.985, 3, 2,1) (0.984, 2, 2,1) (0.961, 3, 3,1) (0.954, 4, 3,1) (0.936, 3, 0, 0) 

(0.929,2,1,0) (0.924,4,4,1) (0.816,2,2,0) (0.677,3,1,0) (0.544,3,2,0) (0.441,3,3,0) 

3 

(0.971,2,2,1) (0.965,3,2,1) (0.961,4,2,1) (0.924,3,3,1) (0.901,4,3,1) (0.851,4,4,1) 

(0.772,2,1,0) (0.620,2,2,0) (0.409,3,-2,0) (0.247,3,-3,0) 

2 

(0.990,2,-1,0) (0.987,2,1,1) (0.981,3,1,1) (0.929,2,2,1) (0.906,4,1,1) 

(0.883,3,2,1) (0.805,3,3,1) (0.765,4,2,1) (0.681,4,3,1) (0.601,4,4,1) (0.150,2,-2,0) 

1 

(0.803,2,1,1) (0.656,2,2,1) (0.148,3,-2,1) (0.088,3,-3,1) 

0.5 

(0.628,2,-2,1) 


TABLE I: Different quadruples yielding the same Qimn- 


additional QNM. 

Nevertheless, some potentially useful bounds may be obtained from detection of a single mode. Suppose we observe 
a Qimn larger than (say) 10. According to Tabled we cannot determine {l,m,n), but we can impose a lower bound 
on j of about 0.88. 

Prospects improve if we assume that we can measure two modes. For definiteness, suppose one mode has Qimn ~ 6 
and the other has Qimn ~ 3. Since they must belong to a quadruple with the same j (within a measurement error of, 
say, one percent), the only possible pairs are, according to Tabled (0.969,4,4,1), (0.929, 2, 2,0), (0.897,4,1, 0) with 
= 6; and (0.971, 2, 2,1), (0.965,3, 2,1), (0.961,4, 2,1), (0.924, 3,3,1), (0.901,4, 3,1) with Qimn = 3. Consider the 
hypothesis that we have detected the pair of modes (0.929, 2, 2,0) with Qimn = 6 and (0.924, 3, 3,1) with Qimn = 3. 
We can test this by computing the mass of the black hole from the two different measured ringdown frequencies. Given 
the measured frequency / 22 o(j = 0.929), we can invert the relation / 22 o(j = 0.929) = 0.703/(27rM) to compute M, and 
then repeat the procedure for the (0.924,3, 3,1) mode. If they yield the same mass, our hypothesis is correct, and the 
measurement is compatible with general relativity. If the masses don’t match, then we may have detected the mode 
(0.969,4,4,1) with Qimn = 6 and any one of the modes (0.971, 2, 2,1), (0.965, 3, 2,1), (0.961,4, 2,1) with Qimn = 3. 
To proceed we take the two different measured values of fimn and we test the compatibility of the resulting masses. 
In this particular example, a determination of the modes involved could be possible. The frequencies corresponding 
to the three different modes, /22i(0.971) ^ 0.79/(27rM), /32i(0.965) ~ 0.94/(27rM), /42i(0.961) ^ 1.13/(27rM), are 
different within the required accuracy. 


V. MULTI-MODE RINGDOWN WAVEFORMS: PRELIMINARIES 

An accurate measurement of QNM frequencies can provide conclusive proof of the astrophysical reality of the black 
hole solutions of general relativity. Most existing studies of ringdown detection 0, S, la 0 S S assume that the 
waveform can be described using only the fundamental (n = 0) I = m = 2 mode of a Kerr black hole. In many 
cases, numerical and perturbative models of astrophysical gravitational wave sources show that, given radiation with 
a certain angidar dependence - that is, given (Z, m) - a good fit of the waveform requires the two lowest modes (n = 0 
and n = 1) p^ . 
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Even more importantly, as stressed by Dreyer et al. 0, a test of the general relativistic no-hair theorem requires 
the identification of at least two QNM frequencies in the ringdown waveform. 

The extension of the formalism to multi-mode situations raises a number of important questions. Which modes 
should we expect to be most relevant in the ringdown waveform? How much energy should we attribute to each QNM 
when a black hole is formed, following either a galaxy merger or the collapse of a supermassive star? Are the energy 
and quality factor of the more rapidly damped modes large enough for them to be detected? Does it make sense to 
talk about the SNR of each QNM, given that (in general) they are not orthonormal in any well-defined sense? Can 
we really discriminate between different QNMs, given that consecutive overtones usually have very similar oscillation 
frequencies? The rest of this paper is devoted to providing preliminary answers to these questions. 


A. Mathematical issues in the definition of mode excitation 

QNMs are relevant to all systems with radiative boundary conditions. For some of these systems QNMs are actually 
a natural extension of an underlying normal mode system. Consider for example the nonradial oscillations of a star. 
The short periods of these oscillations are driven by fluid pressures, and their long damping times are due to the 
(weak) emission of gravitational waves. If we omit gravitational radiation damping we end up with a system that can 
be analyzed in normal modes. In this case we can identify the radiated energy coming from each separate oscillation 
frequency, and decompose the total radiative power into the fraction assigned to each frequency. 

Black hole QNMs (and, for that matter, also the “pure spacetime” w-modes of a star) are different. In this case 
there is a single timescale (given, in geometrical units, by the black hole mass) determining both the frequency and the 
damping time of the oscillations. There is no meaningful way to switch off the radiation damping, and no underlying 
normal mode system. Mathematically, this is reflected in QNMs being eigenfunctions of a non-self-adjoint problem 

QlillilQ. 

This poses difficulties in the definition of a useful and rigorous notion of QNM excitation. In fact, Nollert and Price 
Ellll conjectured that there is no quantitative measure of QNM oscillations satisfying all of the following three 
criteria: (1) the measure is independent of a simple (time) shift of the waveform, (2) the measure can be quantified 
individually for any number of modes, so that the single measures add up to the total norm of the waveform, and (3) 
the measure is useful to quantify the excitation (in particular it lies between, say, 0 % and 100 %). 

Andersson M advocated a more practical viewpoint on this issue. He introduced the following, useful “asymptotic 
approximation”: we require spacetime to be essentially flat in the region of both the observer and the initial data, so 
that initial data should have (compact) support only far from the black hole. Under this assumption we can define 
a mode-decomposition of the time-domain Green’s function which provides an accurate representation of the mode 
excitation and is convergent at late times (see Fig. 2 of and the related discussion). Notice however that we expect 
QNM excitations to arise from data located close to the peak of the perturbative (Zerilli or Regge-Wheeler) potential, 
where spacetime is certainly not flat. For this reason it is not clear how relevant the asymptotic approximation is to 
realistic scenarios. The asymptotic approximation was extended to rotating (Kerr) black holes by Glampedakis and 
Andersson [T^ . Unfortunately they only computed excitation coefficients for scaZor perturbations. 


B. Physical predictions of the energy distribution between different modes 

A rough attempt to deal with a two-mode ringdown waveform in the context of gravitational collapse leading to 
black hole formation can be found in |^. Those authors considered modes with I — m = 2 and Z = 2, m = 0, 
distributing energy between the two through a phenomenological parameter. Is this assumption correct, or should 
other modes be considered as well? How does the energy distribution between modes depend on the physical process 
deforming the black hole? Does the present knowledge of black hole ringdown waveforms provide any information on 
this energy distribution? These are the questions we will try to answer in this Section. 

Ideally, to estimate the relative QNM excitation we would like to have full general relativistic simulations of black 
hole merger and ringdown under different assumptions (different black hole mass ratios, different initial angular 
momenta, realistic initial conditions at the orbital innermost stable orbit). Unfortunately, present state-of-the-art 
simulations in numerical relativity do not provide us with long-term evolutions of black hole mergers, nor with reliable 
estimates of their gravitational wave emission (with the exception of a few, unrealistically symmetric situations). 
Even if we had clean, general relativistic simulations of black hole mergers, there would be additional complications 
of astrophysical nature. Supermassive black hole mergers take place in a very “dirty” galactic environment, and a 
detailed theoretical model of the dynamical interaction between black holes and their surroundings is out of reach, 
given the present understanding of galaxy mergers • 
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Despite these difficulties, we can obtain some insight by considering existing studies of QNM excitation in different 
idealized processes related with black hole formation. The following is a brief summary of interesting results from the 
point of view of black hole excitation. More details can be found in . 


1. Evolution of distorted black holes in full general relativity 

In these simulations, the distortion is sometimes introduced by considering Misner initial data (corresponding to 
two black holes at some given separation), and sometimes evolving initial data corresponding to gravitational waves 
(“Brill waves”). The amplitude of these waves can be large, allowing for the introduction of non-linear effects not 
amenable to perturbation theory. Distorted black hole simulations invariably show that, after a transient depending 
on the details of the initial distortion, quasi-normal ringing dominates the emitted radiation. Most importantly for 
our analysis, they provide some insight into how initial data affect the energy distribution between different QNMs. 
The simulations usually monitor the distortion (more precisely, the ratio of polar and equatorial circumferences) of 
the black hole horizon, fitting the numerical data with the fundamental mode and the first overtone. They typically 
find that the I = 4 horizon distortion is ~ 10“'* smaller than the I = 2 component (see eg. Figs. 2 and 3 in |39l|b 
Quite independently of the initial data and of the black hole spin, the I = 2 component carries away > 95% of the 
gravitational wave energy, but the character of the initial distortion strongly affects the energy distribution between the 
subdominant modes (Table IV of |4f)||b To our knowledge, results from only one nonaxisymmetric (to 0) simulation 
have been produced so far |^. Full three-dimensional simulations of distorted black holes are now computationally 
feasible, and more work in this direction is definitely required. 


2. Simulations of head-on black hole collisions 

Because of its high degree of symmetry, this process has been studied in great detail in full general relativity. 
Existing simulations deal with equal- as well as unequal-mass black holes, either starting from rest or having non¬ 
zero initial momentum (“boosted” black holes). The resulting gravitational waveforms and energy emission are in 
surprisingly good agreement with predictions from linear perturbation theory. These simulations provide at least 
well-motivated and reliable lower bounds on the energy emitted in a realistic merger. 

The present best estimates for equal mass black holes starting from rest predict that the I = 2 component radiates 
~ 0.13% of the black hole mass in gravitational waves [^. Unfortunately wave extraction for the I = 4 component 
has not been carried out yet, and in the past the extraction of this component presented a significant challenge (see 
Fig. 12 of Ref. ^3)- Even if the black holes are not initially at rest, the emitted energy should be less than about 
~ 0.16% of the mass 


3. Black hole formation in gravitational collapse 

With a few exceptions |^ . perturbative and numerical simulations of gravitational collapse usually concentrate 
on stellar mass black holes, so they are not directly relevant to the SMBHs observable by LISA. Nonetheless, some 
predictions of the QNM energy distribution from stellar collapse could carry over (at least qualitatively) to SMBH 
formation induced by the collapse of a supermassive star, and perhaps even to black hole formation following SMBH 
merger. 

Perturbative calculations show that a typical core collapse radiates very little energy (up to ~ 10“^M) in gravi¬ 
tational waves, and that the energy radiated in / = 3 is typically two to three orders of magnitude smaller than the 
I = 2 radiation (see eg. Fig. 9 in |^). 

A classical, nonperturbative axisymmetric simulation by Stark and Piran found that the waveform is very similar 
to the waveform produced by a particle falling into a black hole, and that the radiated energy scales with angular 
momentum j as E/M ~ 1.4x 10“^j"*, saturating at E/M ^ 10“^ for some critical value of j close to the extremal value 
of j = 1. This simulation has recently been extended to the three-dimensional case still keeping a high degree of 
axisymmetry (only m = 0 modes are excited). The new simulations are closer to perturbation theory than the original 
calculation by Stark and Piran, predicting a radiated energy E ~ 1.45 x 10 “®(M/Mq). They also show that the cross 
component of the strain is suppressed by roughly one order of magnitude with respect to the plus component: in our 
notation, they predict that the parameter Nx = Ax/A+ ~ 0.06. We stress again that extrapolation of these results 
to SMBHs is not justified, but the simulations could contain interesting indications on the possible outcome of a full, 
general relativistic black hole merger simulation. 
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VI. MULTI-MODE SIGNAL-TO-NOISE RATIO AND PARAMETER ESTIMATION 

Let us now turn to a multi-mode analysis. For simplicity we assume (as we did in the single-mode case) that spin- 
weighted spheroidal harmonics are real, and we consider a waveform given by the superposition of only two QNMs, 
say 


‘-l-.x 




( 6 . 1 ) 


where we have chosen the phases arbitrarily so that the -I- and x waveforms for both modes are sine functions. We 
denote the mode with indices by a “1”, and the mode with indices {V,m',n') by a “2”. We have seen in 

the single-mode analysis that phases do not play a significant role in parameter estimation, so in principle we could 
simplify things by assuming all phases to be zero. To be able to check and confirm these expectations we only assume 
that = 0, but we allow for a nonzero relative phase of the second mode <() = (t’tm'n' ~ 4’vm'n' Then the 
waveforms take the form 
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Again, the measured waveform is given hy h = h^Fx ■ 

In |22| we show by explicit calculations that, to a good approximation, the angular scalar products between 
spheroidal harmonics for the modes we are interested in are orthonomal on I and m, namely 


Slnin Pi 13'] Si'rn'n'Pi f3]dFl — Si i'dri 


(6.3) 


This suggests that we make a separate treatment of the following two cases: 

(A) Either I p V oi m p m', i.e. we look at modes with different angular dependence. Then the angular scalar 
product between the different modes is zero to a good approximation, and the Fisher matrix can be expressed 
approximately as the sum of the Fisher matrix for mode “1” and the Fisher matrix for mode “2”. This case is 
dealt with in Section Em 

(B) The angular indices I = I' and m = m', but n p n', so we are looking at different overtones with the same 
angular dependence. In this case, as long as we limit attention to the first few modes, the angular scalar product 
is very close to unity for any (n, n') pairs. This case is dealt with in Section IVl Bl 


A. Quasi-orthonormal waveforms (I p I' or m p m') 

We use the FH convention and consider, for simplicity, equal polarization amplitudes Ap = Af = Ai for each mode 
p = 1,2). We further assume the noise to be constant over each mode’s bandwidth. In principle we should assume 
Si p S 21 but if fi and /2 are close enough (which is true in all cases we consider) we can set 5i ~ S '2 — S. Then, a 
simple calculation based on the waveform shows that the total signal-to-noise ratio can be expressed as a sum 
in quadrature. 


2 2,2 

P — Pi + P 2 ) 


(6.4) 


of the single-mode SNRs 

2 _ f \^ Qi 2 _ f MA2 Y Q2piiP P + 2 Ql) 

\ r J 5Tr'^fi{l + 4Ql)S ’ V ^ f2p +‘iQpS ' 

Notice that p 2 becomes symmetric with pi, as it should, in the limit (j) ^ 0. These single-mode SNRs can be obtained 
as appropriate limits of the general expression in the (^-function or constant noise approximation in the FH convention, 
Eq. (13.nil , by choosing the phases there to be cjp = —7r/2, = 0, = —7r/2 -|- (j), and = P- The generalization 
to waveforms involving more quasi-orthonormal modes is straightforward. 

Under our simplifying assumptions, the Fisher matrix for a two-mode waveform depends on five parameters: 
{j, M, Ai, -42, 4>}. For concreteness, we computed the errors on mass and angular momentum obtained by setting 
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FIG. 14: Scaled errors (p<tm, pt^j) for two-mode measurements of mass and angular momentum, with I = m = m' = 2. Solid 
lines refer to l' = 3, dashed lines to l' = 4. 


^ = 0. The analytical expressions are too lengthy to present here. As physical intuition suggests, the crucial element 
in determining the total error is the relative amplitude of the two modes. Without loss of generality, we pick Ai = 1 
and vary the ratio A\/Ai¬ 
ks a second step, we must choose which particular pair of modes we want to compare. We first make the reasonable 
assumption that the dominant modes have n = n' = 0, so we are left with a four-dimensional parameter space, the 
parameters being (/, m, V, m'). Given the standard lore that I = m = 2 modes should in some sense be dominant, we 
examine in detail three cases. 

First, we fix Z = TO = 2, then pick V = m! = i or V = m! = 4. In this case considering additional modes does not 
significantly affect the error with respect to the single-mode case, whatever the amplitude of the second mode. The 
reason is that the functional dependence on j of the mass and angular momentum errors for modes with Z = to is 
basically the same, whatever the value of Z (or to). As we change the relative amplitude of the modes there is a smooth 
transition from, say, the errors on j and M corresponding to Z = to = 2 to the errors corresponding to I' = m' = 3, 
but this transition is almost imperceptible to the eye, the functional behaviors of {p<Jj){j) and {paM/M){j) being 
so similar in the two extreme cases. For this reason we decided not to show any plot, since they would be almost 
indistinguishable from (say) the left panel of Fig. Furthermore, for these corotating modes the fitting functions 
of Appendix ^ do an excellent job at approximating the “true” errors obtained by a numerical calculation of the 
derivatives. 

Secondly, we consider a more physically motivated case in which we fix Z = to = 2, to' = 2, so that both modes 
correspond to a “bar-shaped” deformation. We then look at the effect of setting Z' = 3, 4. Results for these combi¬ 
nations of angular indices are shown in Fig. d From left to right and from top to bottom we assume that the ratio 
of the two mode amplitudes Ai/Ai = {lO^, 10, 1, 10“^, 10“^}. Now the second mode plays some role, and we can 
observe a smooth deformation from the errors corresponding to Z = to = 2 (top left panel) to the errors corresponding 
to Z' = 3, to' = 2 or Z' = 4, m! = 2 (bottom right panel). Notice that the error on mass is roughly independent of Z, 
but, when a single Z dominates, the error on angular momentum scales (roughly) as Z. For Z = to = 2 the mass error 
becomes larger than the angular momentum error at j ~ 0.8. As a consequence of the (rough) scaling with Z of paj, 
this transition moves to larger and larger values of j as the relative amplitude of the second mode grows. When Z' = 3 
dominates, the transition still occurs at j ~ 0.9, but for Z' = 4 the angular momentum error becomes subdominant (if 
at all) only for j > 0.95. 

In Fig. we compute the derivatives and “numerically” (that is, by local interpolating polynomials). 
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FIG. 15: Scaled errors {paM, poj) for two-mode measurements of mass and angular momentum, with I = m = I' = 2. Solid 
lines refer to m! = 1, dashed lines to m! = 0, dot-dashed lines to m! = —1, dotted lines to m' = —2. In the case Ai = 10~^.42 
the error for m' = 0 and m' = — 1 becomes quite large, and for clarity we show results using both a linear and a log scale. 


However, we verified that the fits typically do a very good job, even though they are not as accurate as in the case 
considered above. This is expected, since we are considering corotating modes. For all of these modes does not 
cross zero in the range j € [0,1] (see eg. Fig. EJ, and the fits of Appendix I eI are reasonably accurate. 

The third, and presumably the most physically realistic case, results from fixing I = m = 2, I' = 2 and looking at 
the effect of changing m'. In Fig.^lwe show the smooth deformation of the errors induced by changing Axj Ai in this 
case. Different linestyles correspond to vn! = 1,0, —1, —2, as explained in the caption. The plot is better understood 
as a series of snapshots. The top left panel is an almost-pure I = m = 2 waveform (closely approximating the top 
left panel of Fig. [13 or equivalently the left panel of Fig. nni). Each line becomes more and more dominated by the 
second mode, until (bottom right panels) it approximates quite closely the errors with I = 2 and m = 1,0, —1, —2 we 
displayed in Fig. d In the bottom right panels we clearly see signatures left by the dominance of the second mode. 
For example, the error becomes very large as j —*■ 0 for m! = 0, and as j —> 0.6 for m' = —1. 

In Fig. the derivatives //„„ and Q'l^^ are computed numerically. The fits of Appendix El provide a good 
approximation of the numerical results, except for those cases (i.e. dominance of m' = 0 and m! = —1) where the 
single-mode fits fail to reproduce the location of = 0. 


B. Overtones with the same I and m 


Consider now modes for which I = V and m = m!. To reduce the number of parameters to be estimated, we assume 
that the plus- and cross-components have equal amplitudes (the generalization to unequal amplitudes is trivial). 
Approximating the scalar product between the different spin-weighted harmonics in this case by one, assuming that 
the frequencies are close enough that S*! ~ S '2 — S, and using the expression for the two-mode waveform, we get 
the following SNR: 

2_ 2 , 2 I f A1A2 / 16/i/2Qi(52(/iQ2 +/2Q1) COS0 

P -Pi + P2 + ^-j 



where 


A± = f^Ql + 2/1/2Q1Q2 + Ql [f? + 4(/i ± /2)'Q?] 


(6.7) 
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FIG. 16: Modulation of the total SNR, Eq. 16.61 . induced by the second overtone. We consider an Z = m = 2 perturbation and 
three different values of the angular momentum j (solid lines: j = 0, dashed lines: j = 0.8, dot-dashed lines: j = 0.98). For 
concreteness, we assume that the first overtone has an amplitude A 2 = 10“^Mi. The nearly horizontal (red) lines correspond 
to the SNR in the absence of modulations, Eq. 16.41 . 




FIG. 17: Scaled errors (puM, puj) for the two-mode waveform (solid lines) versus a one-mode waveform (dashed line). Here 
we assume Ai = IOM 2 and consider perturbations with I = m = 2 (left) and Z = 2, m = —2 (right). 


and Pi and p 2 are defined by Eq. 16.511 . 

Fisf.lTi^shows a rescaled SNR given by {r/M)S^^‘^p as a function of (j) for a reasonable amplitude ratio A 2 /A 1 = 10“^. 
The modulation induced by the mixed term is modest, and it scales (roughly) like (.42/-4i) ■ 

In Fig, n 71 we compute the errors on mass and angular momentum for a two-mode waveform with .4i = 10A2j and 
compare the results with a single-mode waveform. We show two representative cases (Z = 2, m = 2 and I = 2, m = —2) 
but we looked at all modes with Z = 2, 3, 4 and all possible values of m. We always found that the correction induced 
by the addition of the second overtone is very small. Indeed, it is so small that it is comparable to the variations 
induced by different prescriptions to compute the error (compare the left panel with the left panel of Fig.lTHll. 

A more realistic treatment would not assume that the second mode has amplitude 1/10 of the first, but would 
attempt to account for the mode excitation using an explicit calculation of the excitations coefficients as functions of 
j. We are currently working on such an approach. 
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VII. RESOLVING QUASI-NORMAL MODES: CRITICAL SIGNAL-TO-NOISE RATIO TO TEST THE 

NO-HAIR THEOREM 


So far we have been discussing multi-mode detection without considering whether a discrimination between two 
different frequencies (in a signal buried in noise) is possible or not. A common rule of thumb to resolve the frequencies 
of two sinusoidal signals with equal amplitudes and quality factors is the so-called Rayleigh criterion (see eg. 
two frequencies /i, /2 are resolvable if 

|/i-/2|At>l, (7.1) 


where At is the duration of the signal. This (purely classical) Rayleigh limit can be beaten, given a signal with 
sufficiently large SNR, as discussed in ^3) ■ Here we introduce a slightly different resolvability criterion that allows 

for different amplitudes {Ai ^ A 2 ) and quality factors (Qi ^ Qn)- In Appendix IDI we will show that actually 

deal with a very special case of our own resolvability criterion. 

To determine whether two quasi-normal mode frequencies are resolvable we first need to determine the measurement 
errors. Here we work in terms of frequency and damping time, instead of M and j. Because Qimn = '^.fimnTimn, the 
Fisher matrix in terms of Timn niay be written, for each index k, = '^fimrJ'kQi^n ^nd T^ = F^ + 

Inverting this Fisher matrix, one can show that the errors in fimn and Timn are given, to leading order 

for large Q, by 


" V27TTP ’ 



(7.2) 


These relations reproduce to leading-order the large-Q behavior predicted by Refs. (which use a different 

parametrization of the waveform), except for a factor in the leading term of cr/. 

Then, for detection of two modes, a natural criterion {d la Rayleigh) to resolve frequencies and damping times is 


1/1 - /2I > max(CT/i,cr/J , |ti - T2\ > max(cr.ri , ) . 


(7.3) 


This means, for example, that the frequencies are (barely) resolvable if “the maximum of the diffraction pattern of 
object 1 is located at the minimum of the diffraction pattern of object 2”. In other words, given two Gaussians with 
variance a we can only distinguish between the peaks if they are separated by a distance larger than a. We can 
introduce two “critical” SNRs required to resolve frequencies and damping times, 


max(pg/,, pcT/J 

1/1-/2I 

max(pcr^^, par2 ) 
\ti-T2\ 

and recast our resolvability conditions as 

p > pcrit = min(p^^i^,pAJ , 
p > pboth = max(p^^;^,pUt). 


Pcrit = 
Pcrit = 


(7.4a) 

(7.4b) 


(7.5a) 

(7.5b) 


The first condition implies resolvability of either the frequency or the damping time, the second implies resolvability 
of both. 

We now need to compute the errors on frequencies (cr/i, cr/ 2 ) and damping times (cTn, <^t 2 ) in a two-mode situation. 
We again use the Fisher matrix formalism. 

We first consider the waveform lES in the quasi-orthonormal case of Sec. IVI AI and for simplicity we also pick 
(/) = 0. Then the analytic expressions for the errors turn out to be very simple: 


P^fi 




1 

2 V 2 


f! (3 + legf) 
AlQl 


' AlQl ' 

./i(l + 4Q2) + /2 (i + 4Q2)_ 



2 

TT 


(3 -I- 4g^) 
AjhQi 


■ AlQl AlQl ' 

./i(l + 4Q2) + /2(i + 4g2)_ 



(7.6a) 


(7.6b) 


Since we consider the “symmetric” case ^ = 0, the errors on /2 and T 2 are simply obtained by exchanging indices 
(1 <-> 2). Notice also that the term in square parentheses is nothing but 57r^S'(r/M)^p^, with p^ given by Eq. (16.411 . 
In the single-mode {A 2 0), large-Q limit, these errors reproduce Eq. (17.211 . 
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FIG. 18: “Critical” SNR pcrit (top) and pboth (bottom) required to resolve different combinations of modes. Here we consider 
the “quasi-orthonormal” perturbations of Sec. IVI Xl leither I ^ I' ov m ^ m'). Following the treatment in that Section, we 
always assume I = m = 2 for the first mode. Different colors/linestyles refer to different values of l' and m' for the second 
mode. Black: I' — 2 and m' = 1 (solid), m' = 0 (dotted), m' = —1, — 2 (dashed, top to bottom). Blue: I' — 3 and m' = 2 
(solid), m' = 3 (dot-dashed). Red: l' = 4 and m' = 2 (solid), m' = 4 (dot-dashed). Different panels refer to different ratios 
A 1 /A 2 , as indicated. 


A calculation of the critical SNRs pcrit and pboth for different QNM pairs is shown in Fig. ^1 Most qualitative 
features of this plot can be explained by looking at Fig.0 First, it is natural that modes with I = I' = 2 are harder 
to resolve than modes with I' ^ 1: in the latter case both frequencies and damping times are different. For the same 
reason, modes with I = V and m ~ m' are harder to resolve. Another predictable feature is that, for modes with 
1 = 1', the Schwarzschild limit is very bad in terms of resolvability: as j —> 0 all frequencies and damping times are 
degenerate with respect to m, so pcrit blows up. In particular, if modes with 1 = 2, m = 2 and 1 = 2, m = 0 are 
dominant (as suggested in [^), resolving them requires large SNRs for black holes with spin j < 0.5, especially if the 
amplitude ratio A 1 /A 2 is not close to unity. The almost-flat j-dependence of pcrit and pboth for modes with m' = 2, 
I' = 3, 4 and V = m! = 3, 4 is in line with the corresponding discussion in Sec. Em the errors on these modes 
have a very similar functional dependence on j, whatever the value of I (or m), so it is natural for the critical SNR 
to be quite insensitive to j. As a rule of thumb, a large rotation parameter j usually helps to resolve modes, mainly 
because of the larger quality factor (remember that we chose the first mode in the pair to have I = m = 2). The only 
exception to this rule is the growth of pboth as j ^ 1 for modes with I = m = 2 and I' = ml(= 3,4). This growth is 
easy to understand: for all modes with I = m the damping tends to zero in the extremal limit, so the denominator of 
Eq. fTlB goes to zero and the corresponding critical SNR blows up. 

We now consider the resolvability of two overtones of modes with the same I and m. Here an explicit calculation of 
the errors on frequency and damping time is not simple, because of the presence of “cross-terms” in the inner products. 
The final expressions are much more lengthy and involved than Eq. so we do not include them here. Instead 

we carried out the calculation of Eqs. - (I7.5bll numerically. We assumed that the waveform is a superposition of 

the fundamental mode and the first overtone with some given values of (I, m), and made the plausible educated guess 
that the amplitude of the overtone is a factor ten smaller than the amplitude of the fundamental mode {Ai = 10 ^ 2 ), 
meaning that the energy carried by the first overtone is roughly 1% of the total energy radiated in the ringdown. 

Our results for pcrit and pboth are shown in Figs.|niand01 respectively. Once again, most features of these plots can 
be explained by looking at Fig.[3 As long as we are content with resolving either the frequency or the damping time, 
higher Z’s and positive m’s are easier to resolve because the frequency separations are similar, but the quality factor is 
larger. In fact, as j grows the quality factor for corotating modes grows, but it decreases for counterrotating modes. 
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and this affects resolvability in the corresponding way. If we want to distinguish both frequency and damping time, 
the behavior of the critical SNR changes. The corotating fundamental mode with I = m is impossible to resolve from 
the first overtone in the extremal limit j 1 because they become degenerate, tending to the Detweiler frequency. 

Fig.El shows that, quite independently of (I, m) and of the angular momentum j, resolving either the frequency 
or the damping time requires a SNR p > pcrit 10^. Since the addition of a small-amplitude overtone does not 
significantly alter the SNR, we can use the SNR predictions of Fig. 13 to deduce that tests of the no-hair theorem 
should be feasible even under the most pessimistic assumptions on the ringdown efficiency Crd (at least for equal-mass 
mergers) as long as the first overtone radiates a fraction ~ 10“^ of the total ringdown energy. However, resolving 
both frequencies and damping times typically requires a SNR p > pboth ~ 10^. This is only possible under rather 
optimistic assumptions on the radiative efficiency, and it can be impossible if the dominant mode has I = m = 2 and 
the black hole is rapidly spinning (solid black line in the left panel of Fig.^. We hope to refine this analysis using 
computed estimates of excitation coefficients. 


VIII. CONCLUSIONS 

We have presented a general framework for analysing the detectability of quasi-normal ringdown gravitational waves 
from massive black holes, and for using them to estimate parameters of the hole and to test the general relativistic 
no-hair theorem. In this initial work, we made a number of simplifying assumptions, including the reality of the 
spin-weighted spheroidal harmonics (thus simplifying the angle averaging), the restriction to one-mode or two-mode 
situations, and the assumption of large quality factors in deriving analytic expressions. In future work, we hope to 
explore the consequences of relaxing some of these assumptions. Although we focused on detection of modes using 
LISA, our framework is equally applicable to ground-based interferometers. 

For example, the non-angle-averaged case can be studied using Monte Carlo simulations in place of angle averaging, 
thereby permitting inclusion of the fully complex spheroidal harmonics. Multi-mode calculations involving more than 
two modes should be carried out. Our discussion of testing the no-hair theorem has focused only on resolvability 
of modes; we hope to use our tools to perform in detail the no-hair test suggested in Q. A crucial point is to 
determine the excitation of modes. When they have different angular dependence, we must use numerical relativistic 
simulations of mergers, or minimally simulations of distorted black holes with initial data that mimics the merger. 
For overtones with the same angular dependence, we can use perturbative and analytical techniques to estimate the 
relative amplitudes; this work is in progress. 

It will be important to extend the analysis from matched-filtering to more sophisticated data-analysis techniques 
(eg. the tiling method used by the TAMA group @,^3) combine them with time-delay interferometry (TDI), 
especially in the case of short damping times. 

Another important question is the effect of combining parameter estimation from ringdown with parameter esti¬ 
mation from inspiral. This may make it possible to improve the measurement of mass or angular momentum by 
taking advantage of “prior” information, and may permit one to deduce how much energy and angular momentum is 
radiated in the merger phase. 
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APPENDIX A: ACCURACY OF THE 5-FUNCTION APPROXIMATION 

How accurate is the i5-function approximation we employed in Sec. nnn' The answer to this question depends on 
{l,m,n) and on the angular momentum j of the black hole. In general, we expect the approximation to hold better 
in the high-Qj^n limit: that is, for low overtones and (when the black hole is spinning) for modes with m > 0. This 
expectation is made more quantitative in Fig. 1191 where we plot the relative error on the SNR due to the (5-function 
approximation for different modes as a function of the black hole mass, for D/, = 3 Gpc. To compute the “exact” SNR, 
we use the FH SNR in the Breit-Wigner form [Eq. 13.911 or I3.12all ]. The plot refers to a near-extremal black hole 
(j = 0.98), but the qualitative features we observe would be the same for different (non-zero) values of the rotation 
rate. In this near-extremal case the numerical integration of the “full” SNR becomes problematic for corotating modes 
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FIG. 19: Modulus of the error in the SNR due to the d-function approximation, \{ps — pexact)/pexact|, for modes with I = 2 
and I = 4, for j = 0.98 and Dl = 3 Gpc. For each I we consider the fundamental mode (n = 0) and the second overtone 
(n = 2). The horizontal line marks a 10% deviation from the “exact” result: for those modes lying below the horizontal line the 
5-function approximation can be considered very accurate. The “inverted peaks” are just zero-crossings of |(p 5 — pexact)/pexact|. 


(to > 0), because the Breit-Wigner function effectively approaches a (^-function. This problem usually shows up for 
very large masses (this is the origin of the small wiggles in the high-mass tail of the curve marked by to = 2 in the 
top-left panel). To achieve some required numerical accuracy it is sufficient to double the number of integration points 
until the high-mass wiggles disappear and the integration routine converges. 

The plot confirms our expectations. As a general trend, the (5-function approximation is more accurate for corotating 
modes (to > 0), for the simple reason that the quality factor Qimn of corotating modes tends to infinity (or becomes 
very large) as j —> 1. The approximation is more accurate for slowly damped modes (small values of the overtone 
index n), which have a larger quality factor. As a rule of thumb, we can say that the (5-function approximation is 
reasonably accurate whenever it induces deviations smaller than ^ 10% with respect to the exact result: all modes 
lying below the horizontal lines in Fig. 1191 satisfy this criterion. Fig. llOl shows that the approximation is very accurate 
for the fundamental (n = 0) mode, whatever the values of (l,m). However, when we consider the n = 2 overtone, it 
is only marginally accurate for corotating modes, and inaccurate for modes with to < 0. 

The large deviation at mass M ^ 10^Mq corresponds to the frequency band where white-dwarf confusion noise 
dominates over instrumental noise. In this frequency band the noise curve is not a very smooth function, and the 
SNR (especially for corotating modes) is more sensitive to the “full” shape of the noise curve. We verified that 
the peak disappears when we omit the white-dwarf confusion noise from the Barack-Cutler noise curve. The larger 
disagreement in this frequency regime corresponds, unfortunately, to the most promising mass range for detection 
of black hole ringdown from galactic centers. This means that we probably need a good control of the white-dwarf 
confusion noise (and a rather accurate model of the ringdown waveform) to detect these events. 
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APPENDIX B: FISHER MATRICES FOR A SINGLE-MODE WAVEFORM 

Taking the Fourier transform of ringdown waveforms requires some care. For ringdown waves from SMBH mergers, 
the waveform emitted before ringdown starts (say, for t < 0) is the poorly known merger waveform. To obtain the 
SNR from the ringdown signal alone, a possible guess is to assume that the waveform h{t) = 0 for t < 0. Under this 
assumption we can calculate the Fourier transform of the waveform by integrating only over the range t > 0. This 
method was used, among others, by Echeverria and Finn |^, and for this reason we call it the Echeverria-Finn 
(EF) convention. More rigorously, we would really like to know the probability that a ringdown waveform is present in 
the data stream starting (say) at t = 0. This probability should be computed integrating over all possible realizations 
of the noise for t < 0. Appendix A of FH shows that this is equivalent to minimizing the SNR over all choices of the 
function h(t) on the negative t axis, and that the SNR obtained by minimizing over these choices is always within a 
few tens of a percent of the SNR obtained using the “FH prescription”, as described in Section llTT A I In the body of 
the paper we usually adopt the FH prescription. In this Appendix we present, for completeness, the Fisher matrix 
calculation in the EF convention. We also list the elements of the Fisher matrix adopting the FH convention, but 
assuming that we do not know a priori the relative amplitudes or phases of the two polarizations. 


1. Echeverria-Finn convention, four parameters 


Here we compute the single-mode Fisher matrix following the Echeverria-Finn convention for calculating the Fourier 
transforms. We start with the waveform as defined in Eqs. (lOll - (lOll . But in this case, instead of taking \t\ in the 
damped exponential, integrating over t from — oo to -boo and dividing the resulting spectrum by -\/2 to compensate 
for the doubling, we assume that the waveform vanishes for t < 0 , and integrate only over positive t. 

Assuming that the noise can be considered constant over the bandwidth of the signal, we get the angle-averaged 
SNR pef = 7(1 + 'iQfmn ~ P + ‘^Qimnd), where a, (3 and 7 are given by Eqs. (14.1011 . The result is equivalent to 
Eq. E2S1). 

The Eisher matrix in the (A+, fimn, Qimn) basis is given by 



— ^^- 1-^2 (1 4” 4(5imn P3-2QimnC() , 

(Bla) 


= (Q^ + 2QlmnP) , 

(Bib) 


“ 2 A+ fi ^Qlmn P 4" ‘^QlmnCp) j 

(Blc) 

^A+Qi^^ 

— 94 + n l_l_zLn 2 ‘^Qlmn) (1 '^QlmrPP^Ql^nno] , 

(Bid) 

^<l>Ln4>Lr. 

= 7 ( 14 - ^Qlmn 3- P — ‘2QlmnO:) , 

(Ble) 

F ,+ f 

= -^7^[«-2Qimn(l+4QLn-/3)] , 

^Jlmn 

(Blf) 


— 9 ^ 1 4Q;„„)a-b 4Qimn/3] , 

^Wlmn \^^^'^lmn/ 

(Big) 


r/imn/imn “ 9^2 /3 + 2QlmnOi\ , (Blh) 

P Imn 


^Qlmr^Ql 


1 _ 

Vi mn Qi mn 1 -|- 4Q 

7 _ 1 

2 QLn(i + 4 gLJ" 


2 - [(1 + 

Imn 

[(1 + ‘3Q}mnf ~ 


2 - (1 - 4QLn)/3 + 4Qz„„a] , 

(1 — 12Qi^j^)P -|- (6 — 8Qi^j^)QlmnCx] ■ 


(Bli) 

(Blj) 













33 


In the (A+, M, j) basis, an expansion of the errors in powers of gives 
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For the correlation coefficients we get 
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2. Flanagan-Hughes convention, six parameters 

Here we list the elements of the Fisher matrix for the case where we do not know a priori the relative amplitudes 
or phases of the two polarizations. In the six-parameter basis (H+, 40rnn^ fimn, Qimn), they are 
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= 94 f (1 + 4gL« + cos2^/^„) cost/-, 

(B4f) 

^flmn 

= 94 f (l + 4gL« cos2(j)0^0 sm0 , 

Imn 

(B4g) 

A'^ Qiyuyi 

= 94 n 1 -C 4(02 [(l+4gLn)^ + (l 4QL«) cos 2</)+„„] cost/-, 

2AQimn 1 + ^Qljnn 

(B4h) 

^A^ Qlmn 

= 94n 1-C4n2 [(l + 4gLn)^ (1 4QLn)cOs2()/fJ SIUt/-, 

^^Wlran t -|- 

(B4i) 


= 7 (I -f 4Q^^„ - COS 2(t0) cos^ t/- , 

(B4j) 

r < X -X 

= 7 (I -b 4gL„ -b cos 2<l)^ ) sin^ t/- , 

(B4k) 

^ Imn^ Iran 

= 0, 

(B41) 
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r/,™„/,„. 

^flmr^Qlm^ 

^Qlmr^Qlm 


7 


Vi mn 

7 


2flr. 


COS^V', 

-sin2(/)f„„ sin2^/>, 


7 


1 ^Qfmn 

1 + iQ 


2Qi mn 

7 fi-^Ql 


Imn 


2Qi mn 
7 

‘^flmn 


1 + 4 QLn 


sin 2(j}'^ cos^ '0 ) 


sin 20 ^ sin^ ip , 


[(1 + — /?] , 

7 1 


2 /i mn Qi mn 1 + ‘^Qimn 
7 1 


[(l + 4QLnf-(l-4QLJ/3] , 


2Q2 (1 + 4Q2 )2 

^ Imn V ' ^ Imny 


[(l + 4gLnf-(1-12QLJ/3] , 


(B4m) 

(B4n) 

(B4o) 

(B4p) 

(B4q) 

(B4r) 

(B4s) 


where = (A+)^ + cos0 = AA jA and sin0 = A^ jA, and where a, /3 and 7 are given by Eqs. (14.1011 . Also, 

as in the four-parameter case, p]pjj = 7(1 + 4(5^^„ — j3). 

As before, one can convert to the {M,j) basis using Eqs. 114.1111 . and then invert the Fisher matrix to obtain the 
errors and correlation coefficients. The resulting formulae are lengthy and unenlightening, and we do not display them 
here. 


APPENDIX C: LISA NOISE CURVE 

Generally, the LISA community has been using the so-called sky-averaged spectral noise density 5®^ [see e.g.. 
Ref. and the LISA Pre-Phase A Report], computed by a combination of three factors, including: (i) the raw spec¬ 
tral noise density Sn , (ii) the gravitational-wave transfer (response) function R and (iii) the noise transfer (response) 
function i?„. They combine together in |5l| 


CySA _ ^ri Rn 


(Cl) 


In the low frequency limit, the GW transfer function used in the LISA Sensitivity Curve Generator is i? = 
4('\/3/2)^l/5 = 3/5, where the factor (-03/2)^ comes from the LISA arms being at 60°, the factor 1/5 is due to the 
sky-average of the pattern functions = 1/5) and the factor 4 depends on the particular read-out variable used. 

Since our definition of the GW signal already includes the factor •\/3/2, and since we choose to sky-average the signal, 
we use an effective non-sky-averaged spectral density, obtained by multplying by (-\/3/2)^/5 = 3/20. The final 
result is Eq. 113. 2011 ■ and has been obtained also in Ref. 0|- We estimate white-dwarf confusion noise following Qj, 
which uses results from l5^ : the galactic contribution is approximated as 


(/) = 2.1 X 10-45 ' Hz-4, 


(C2) 


and the contribution from extra-galactic white dwarfs as 

^ex-gal{^) ^ 4 2 X 10-47 ' Hz-^ . 


(C3) 


We compute the total (instrumental plus confusion) noise as 

5^(/) = min {5r^(/)/exp {-I df) , S^^if) + '(/)} + ■ 


(C4) 


Here dN/df is the number density of galactic white-dwarf binaries per unit gravitational-wave frequency, for which 
we adopt the estimate 


— = 2 X 10-5 Hz 

df 


,-i 


1 Hz y 4/3 


/ 


(G5) 
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FIG. 20: Analytic approximation to the LISA root noise spectral density curve used in this paper and in Ref. 0 (dashed line) 
and the curve produced using the LISA Sensitivity Curve Generator (solid line). The SCG curve has been multiplied by 
a factor of -^3/20 to obtain an effective non-sky averaged noise spectral density. The SCG noise curve does not include the 
extragalactic white dwarf confusion noise while the analytical approximation curve does. 


“ ^mission discretely Fourier transformed data for a LISA mission lasting a time Tmission; 

and K ~ 4.5 is the average number of frequency bins that are lost when each galactic binary is fitted out. The factor 
exp (—thus represents the fraction of “uncorrupted” bins where instrumental noise still dominates. 
The analytic root noise spectral density curve El used in this paper is shown in Fi g. |2()l together with the corre¬ 
sponding root noise spectral density curve from the LISA Sensitivity Curve Generator |^. The SCG curve shown is 
obtained using the nominal values SNR=1, arm length = 5 x 10® m, telescope diameter = 0.3 m, laser wavelength 
= 1064 nanometers, laser power = 1.0 Watts, optical train efficiency = 0.3, acceleration noise = 3x 10“^^ ms“® 
and position noise budget = 2 x 10~^^ mHz ~^^®, with position noise setting the floor at high frequency. The data 
returned by the SCG is then multiplied by y/3/20 to obtain the effective non-sky averaged curve shown in Fig. [201 


APPENDIX D: THE SHAHRAM-MILANFAR CRITERION FOR RESOLVABILITY 


For damped sinusoids it makes sense to assume that the signal duration At is given by the e-folding time r, so the 
“standard” Rayleigh criterion GU becomes (/i —/ 2 )t > 1. In 0 ^1®° E3) authors discuss how to go beyond 

this (purely classical) Rayleigh limit. They show how to bypass this limit for two sinusoids with nearby frequencies 
fi ) /2 such that 2(5 = /i — /2 ^/B, where B is the duration of the signal. Assuming white noise with zero mean 

and unit variance cr®, th ey g et a condition for the minimum SNR required for the frequencies of two equal-amplitude 
sinusoids to be resolved 1^: 


640 {T-\Pf)-T-\Pa)f 
7r4 {25BY 


(Dl) 


Here T{x) is the right-tail probability function for a Gaussian random variable X with zero mean and unit variance 
[ 5 ^. Pd and Pf denote the detection and false-alarm rates, typical values for these being, say, Pd = 0.99 , Pf = 0.01. 
Their analysis can be carried over to our case, assuming large Q factors, and therefore a signal with duration r ~ Q/ f- 
We would get 


640 (r-i(P/)-T-i(Pd))' 1 

r\h-hY ^t4(/i-/2)4- 


(D2) 


For overtones (“nearly-parallel signals”, in the sense of the scalar product of SWSHs being close to unity) with 
equal amplitudes and quality factors the frequency error has a large-Q behavior 

1 


- / 2 I 


(D3) 
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TABLE II: First three overtones for I — 2. A comma separates the real part from the imaginary part of Mui. To save space, 
in this and the following Tables we omit leading zeros. 

I = 2, n = 0 

j m = 2 m = l m — 0 m = —1 m = —2 
0.00 .3737,.0890 .3737,.0890 .3737,.0890 .3737,.0890 .3737,.0890 
0.10 .3870,.0887 .3804,.0888 .3740,.0889 .3678,.0890 .3618,.0891 
0.20 .4021,.0883 .3882,.0885 .3751,.0887 .3627,.0889 .3511,.0892 
0.30 .4195,.0877 .3973,.0880 .3770,.0884 .3584,.0888 .3413,.0892 
0.40 .4398,.0869 .4080,.0873 .3797,.0878 .3546,.0885 .3325,.0891 
0.50 .4641,.0856 .4206,.0862 .3833,.0871 .3515,.0881 .3243,.0890 
0.60 .4940,.0838 .4360,.0846 .3881,.0860 .3489,.0876 .3168,.0890 
0.70 .5326,.0808 .4551,.0821 .3941,.0845 .3469,.0869 .3098,.0887 
0.80 .5860,.0756 .4802,.0780 .4019,.0822 .3454,.0860 .3033,.0885 
0.90 .6716,.0649 .5163,.0698 .4120,.0785 .3444,.0849 .2972,.0883 
0.98 .8254,.0386 .5642,.0516 .4223,.0735 .3439,.0837 .2927-.0881 


_ 1^2, n = l _ 

j m = 2 m = l m — 0 m = —1 
0.00 .3467,.2739 .3467,.2739 .3467,.2739 .3467,.2739 
0.10 .3619,.2725 .3545,.2731 .3472,.2737 .3400,.2744 
0.20 .3790,.2705 .3635,.2717 .3486,.2730 .3344,.2744 
0.30 .3984,.2680 .3740,.2698 .3511,.2718 .3296,.2741 
0.40 .4208,.2647 .3863,.2670 .3547,.2700 .3256,.2734 
0.50 .4474,.2602 .4009,.2631 .3594,.2674 .3225,.2723 
0.60 .4798,.2538 .4183,.2575 .3655,.2638 .3201,.2708 
0.70 .5212,.2442 .4399,.2492 .3732,.2585 .3184,.2686 
0.80 .5779,.2281 .4676,.2358 .3826,.2507 .3173,.2658 
0.90 .6677,.1953 .5059,.2097 .3935,.2385 .3167,.2620 
0.98 .8249,.1159 .5477,.1509 .4014,.2231 .3164,.2581 


m = —2 
.3467,.2739 
.3330,.2750 
.3206,.2759 
.3093,.2765 
.2989,.2769 
.2893,.2772 
.2803,.2773 
.2720,.2773 
.2643,.2772 
.2570,.2770 
.2515,.2768 


j m = 2 
0.00 .3011,.4783 
0.10 .3192,.4735 
0.20 .3393,.4679 
0.30 .3619,.4613 
0.40 .3878,.4533 
0.50 .4179,.4433 
0.60 .4542,.4303 
0.70 .4999,.4123 
0.80 .5622,.3839 
0.90 .6598,.3275 
0.98 .8238,.1933 


_ 1 = 2, n = 2 _ 

m = l m = 0 m = —1 m = —2 
.3011,.4783 .3011,.4783 .3011,.4783 .3011,.4783 
.3104,.4756 .3017,.4778 .2932,.4801 .2846,.4825 
.3214,.4719 .3038,.4764 .2866,.4811 .2697,.4862 
.3342,.4671 .3074,.4739 .2813,.4814 .2559,.4895 
.3492,.4607 .3124,.4701 .2772,.4808 .2433,.4925 
.3669,.4522 .3190,.4647 .2741,.4794 .2316,.4952 
.3878,.4407 .3273,.4571 .2721,.4768 .2207,.4977 
.4133,.4241 .3374,.4464 .2709,.4729 .2107,.4999 
.4451,.3984 .3488,.4307 .2703,.4674 .2013,.5019 
.4867,.3502 .3591,.4067 .2697,.4600 .1925,.5038 
.5201,.2331 .3599,.3808 .2686,.4527 .1858,.5051 


so our frequency resolvability criterion, |/i — fpl > max (cr/j, cry^) agrees with the Shahram-Milanfar result (ID2II . 
Notice however that the situation considered in |48| is very special: overtones with the same equal amplitudes 

A\ = A 2 and equal quality factors Qi = Q 2 are the only class of modes for which l/|/i — / 2 |- 


APPENDIX E: QUASINORMAL FREQUENCIES FOR ROTATING BLACK HOLES 

There is a vast literature on Kerr QNMs IH, l6fi| , but the QNM frequencies which are relevant for detection 
have never been systematically tabulated. Here we list for reference the values of the complex frequencies and 
separation constants for selected QNMs. These were calculated using two different numerical implementations of the 
continued fraction method of Leaver |57|| . one written by one of us (EB), and the other provided by Hisashi Onozawa. 
Both methods give results in excellent agreement with each other. Shown in Tables im 11111 and HU are the real and 
imaginary parts of the complex QNM frequencies, separated by commas, for Z = 2, 3 and 4, respectively, and for 
the three lowest overtones. Tables Ivl IVII and lyill display the angular separation constants Aimn in the same format. 
Following Leaver, the numerical values of the Aimn’s were obtained assuming that the perturbations have a time 
dependence of the form while in this paper we use the opposite convention on the Fourier transform. 

The numerical values of the quasi-normal frequencies we listed in Tables HU EH and Hvl can be fitted to reasonable 
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TABLE III: First three overtones for Z = 3. 





0 

II 

e 

co" 

II 




J 

CO 

II 

m = 2 

m = 1 m = 0 

m = —1 

m = —2 

m = —3 


0.00 .5994,.0927 .5994,.0927 .5994,.0927 .5994,.0927 .5994,.0927 .5994,.0927 .5994,.0927 

0.10 .6208,.0924 .6137,.0925 .6067,.0926 .5999,.0926 .5932,.0927 .5867,.0928 .5802,.0928 

0.20 .6448,.0920 .6297,.0921 .6153,.0923 .6014,.0924 .5880,.0926 .5752,.0927 .5628,.0929 

0.30 .6721,.0913 .6480,.0915 .6252,.0918 .6038,.0921 .5837,.0923 .5647,.0926 .5469,.0928 

0.40 .7037,.0902 .6689,.0906 .6369,.0911 .6074,.0915 .5802,.0920 .5553,.0924 .5323,.0927 

0.50 .7409,.0887 .6934,.0893 .6506,.0900 .6121,.0908 .5776,.0915 .5467,.0920 .5188,.0925 

0.60 .7862,.0864 .7228,.0873 .6670,.0885 .6183,.0897 .5758,.0908 .5388,.0916 .5063,.0922 

0.70 .8437,.0829 .7592,.0842 .6870,.0861 .6261,.0882 .5749,.0899 .5316,.0912 .4946,.0919 

0.80 .9219,.0770 .8068,.0791 .7121,.0824 .6360,.0860 .5749,.0888 .5250,.0906 .4837,.0916 

0.90 1.0446,.0655 .8762,.0689 .7455,.0758 .6487,.0826 .5758,.0872 .5191,.0900 .4735,.0913 

0.98 1.2602,.0387 .9769,.0453 .7833,.0643 .6615,.0782 .5773,.0856 .5146,.0894 .4657,.0910 


Z = 3, n = 1 


J 


m = 3 


m — 2 


m = 1 


m = 0 


m = —1 


m = —2 


m = —3 


0.00 .5826,.2813 .5826,.2813 .5826,.2813 .5826,.2813 .5826,.2813 .5826,.2813 .5826,.2813 

0.10 .6053,.2802 .5978,.2805 .5904,.2808 .5832,.2811 .5761,.2814 .5691,.2817 .5623,.2820 

0.20 .6306,.2785 .6148,.2791 .5995,.2798 .5848,.2804 .5706,.2811 .5569,.2817 .5437,.2823 

0.30 .6593,.2761 .6341,.2771 .6102,.2782 .5876,.2793 .5662,.2804 .5459,.2814 .5268,.2823 

0.40 .6924,.2726 .6563,.2741 .6227,.2758 .5916,.2776 .5627,.2793 .5360,.2808 .5112,.2822 

0.50 .7312,.2676 .6821,.2698 .6375,.2724 .5969,.2751 .5602,.2777 .5270,.2800 .4967,.2818 

0.60 .7782,.2604 .7130,.2634 .6550,.2674 .6038,.2717 .5587,.2756 .5188,.2789 .4833,.2813 

0.70 .8375,.2494 .7510,.2538 .6762,.2600 .6124,.2668 .5580,.2728 .5113,.2775 .4708,.2807 

0.80 .9176,.2315 .8004,.2379 .7025,.2484 .6230,.2598 .5582,.2692 .5045,.2758 .4590,.2800 

0.90 1.0425,.1966 .8714,.2068 .7362,.2277 .6357,.2491 .5593,.2644 .4983,.2739 .4480,.2792 

0.98 1.2599,.1161 .9706,.1341 .7679,.1930 .6468,.2359 .5606,.2593 .4938,.2720 .4396,.2784 






1 = 3, n 

= 2 




j 

m = 3 

m = 2 

m = 1 m 

= 0 

m = —1 

m = —2 

m = —3 


0.00 .5517,.4791 .5517,.4791 .5517,.4791 .5517,.4791 .5517,.4791 .5517,.4791 .5517,.4791 

0.10 .5766,.4763 .5684,.4771 .5603,.4779 .5523,.4787 .5445,.4795 .5368,.4803 .5292,.4812 

0.20 .6043,.4725 .5872,.4741 .5705,.4758 .5543,.4775 .5386,.4792 .5234,.4809 .5086,.4826 

0.30 .6356,.4674 .6084,.4699 .5825,.4725 .5577,.4753 .5340,.4781 .5114,.4809 .4897,.4836 

0.40 .6714,.4605 .6328,.4640 .5966,.4679 .5625,.4721 .5305,.4763 .5005,.4804 .4723,.4842 

0.50 .7131,.4511 .6612,.4558 .6131,.4614 .5689,.4675 .5282,.4736 .4907,.4793 .4561,.4845 

0.60 .7632,.4379 .6948,.4441 .6327,.4522 .5769,.4611 .5270,.4698 .4819,.4778 .4410,.4846 

0.70 .8258,.4185 .7358,.4268 .6561,.4387 .5868,.4521 .5267,.4648 .4739,.4757 .4269,.4843 

0.80 .9096,.3874 .7884,.3988 .6844,.4177 .5984,.4391 .5272,.4581 .4667,.4730 .4136,.4839 

0.90 1.0383,.3282 .8622,.3449 .7177,.3806 .6106,.4197 .5282,.4492 .4600,.4696 .4011,.4834 

0.98 1.2592,.1935 .9605,.2181 .7349,.3224 .6176,.3976 .5286,.4403 .4550,.4665 .3916,.4828 


accuracy by simple functions, which can easily be inverted to yield j , M once the ringing frequencies are known. The 
fitting functions all have the form 


^iTTin = MuJimn = /l + / 2(1 “ , 

Qlmn = gi+g2(l-j)®- 


(El) 

(E2) 


In Tables IVInlIIXI and IXI we list the fitting coefficients fi (for the frequency) and qi (for the quality factor). 

Our fits (and the quoted values for the errors) refer to the range j G [0,0.99]. The only such fits carried out h 
other authors to date treated only the I = 2, m = 2 and I = 2, m = 0 modes. For I = 2, m = 2 Echeverria 
proposed the fit 


^220 - [1-0.63(1-j)°'3] , 
Q 220 ^ 2 (l-j)- 0 - 45 . 

For I = 2, m = 0 Fryer, Holz and Hughes found 


(E3) 

(E4) 


^200 - - [1 - 0.13(1-j)° ®] 
Q 200 — 3 — (1 — j)°'^ . 


(E5) 

(E 6 ) 
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j m = 4 
0.00 .8092,.0942 
0.10 .8387,.0940 
0.20 .8717,.0935 
0.30 .9092,.0929 
0.40 .9525,.0918 
0.50 1.0034,.0903 
0.60 1.0650,.0879 
0.70 1.1427,.0842 
0.80 1.2475,.0781 
0.90 1.4104,.0662 
0.98 1.6919,.0388 


m = 3 
.8092,.0942 
.8313,.0940 
.8560,.0936 
.8839,.0930 
.9159,.0921 
.9532,.0907 
.9978,.0885 
1.0530,.0852 

1.1254.. 0796 

1.2322.. 0685 

1.3944.. 0425 


j m = 4 
0.00 .7966,.2843 
0.10 .8271,.2835 
0.20 .8611,.2821 
0.30 .8997,.2799 
0.40 .9441,.2766 
0.50 .9962,.2717 
0.60 1.0591,.2645 
0.70 1.1382,.2532 
0.80 1.2445,.2347 
0.90 1.4090,.1986 
0.98 1.6917,.1165 


m = 3 
.7966,.2843 
.8194,.2837 
.8449,.2824 
.8738,.2805 
.9067,.2775 
.9451,.2731 
.9909,.2665 
1.0474,.2562 

1.1212.. 2391 

1.2293.. 2055 

1.3917.. 1268 


j 

m = 4 

m = 3 

0.00 

.7727,.4799 

.7727,.4799 

0.10 

.8049,.4780 

.7969,.4784 

0.20 

.8409,.4751 

.8239,.4759 

0.30 

.8815,.4708 

.8544,.4721 

0.40 

.9280,.4647 

.8892,.4666 

0.50 

.9824,.4559 

.9296,.4586 

0.60 

1.0478,.4431 

.9775,.4469 

0.70 

1.1295,.4237 

1.0365,.4291 

0.80 

1.2387,.3921 

1.1129,.3998 

0.90 

1.4061,.3313 

1.2236,.3428 

0.98 

1.6913,.1942 

1.3871,.2096 


TABLE IV: First three overtones for I = 4. 


i = 4, n = 0 


m = 2 


m = 1 


m = 0 


m = —1 


m = —2 


m = —3 m = —4 


.8092,.0941 
.8240,.0940 
.8407,.0937 
.8597,.0932 
.8815,.0924 
.9069,.0912 
.9369,.0893 
.9734,.0865 
1.0198,.0820 
1.0836,.0733 
1.1631,.0558 


.8092,.0942 
.8168,.0941 
.8259,.0938 
.8366,.0934 
.8493,.0927 
.8642,.0917 
.8820,.0902 
.9035,.0881 
.9301,.0848 
.9647,.0791 
1.0019,.0703 


.8092, 

.8098, 

.8116, 

.8146, 

.8191, 

.8250, 

.8326, 

.8423, 

.8546, 

.8702, 

.8860. 


.0942 

.0941 

.0939 

.0935 

.0930 

.0922 

.0911 

.0896 

.0874 

.0841 

.0800 


.8092, 

.8028, 

.7977, 

.7937, 

.7908, 

.7890, 

.7883, 

.7889, 

.7907, 

.7941, 

.7979, 


0942 

0941 

0940 

0937 

0933 

0927 

0919 

0909 

0895 

0877 

0858 


.8092, 

.7960, 

.7842, 

.7737, 

.7643, 

.7559, 

.7485, 

.7420, 

.7363, 

.7314, 

.7282, 


0942 

0942 

0941 

0938 

0935 

0931 

0925 

0919 

0911 

0901 

0892 


.8092, 

.7893, 

.7712, 

.7546, 

.7395, 

.7255, 

.7125, 

.7005, 

.6894, 

.6790, 

.6712, 


0942 

0942 

0941 

0940 

0937 

0934 

0930 

0926 

0921 

0915 

0910 


.8092, 

.7826, 

.7585, 

.7365, 

.7162, 

.6974, 

.6800, 

.6637, 

.6485, 

.6343, 

.6235, 


.0942 

.0942 

.0942 

0941 

0939 

0936 

.0933 

.0930 

0926 

.0922 

0919 


m = 2 
.7966,.2843 
.8119,.2838 
.8293,.2828 
.8489,.2811 
.8715,.2785 
.8977,.2747 
.9287,.2691 
.9663,.2604 
1.0137,.2464 
1.0779,.2199 

1.1519,.1663 


I — 4, n — 1 


m = 1 
.7966,.2843 
.8046,.2840 
.8140,.2831 
.8252,.2817 
.8384,.2796 
.8539,.2765 
.8725,.2720 
.8947,.2653 
.9221,.2551 
.9565,.2376 
.9899,.2113 


m = 0 
.7966,.2843 
.7973,.2841 
.7992,.2835 
.8025,.2823 
.8072,.2806 
.8136,.2782 
.8217,.2748 
.8319,.2701 
.8446,.2634 
.8602,.2532 
.8746,.2407 


m = —1 
.7966,.2843 
.7901,.2843 
.7848,.2838 
.7808,.2829 
.7780,.2816 
.7763,. 2799 
.7759,.2775 
.7767,.2743 
.7788,. 2702 
.7822,.2646 
.7857,.2586 


m = —2 
.7966,.2843 
.7830,. 2844 
.7709,.2841 
.7601,.2835 
.7505,.2826 
.7420,.2813 
.7345,.2797 
.7279,.2777 
.7222,.2752 
.7174,.2723 
.7140,.2694 


m = —3 
.7966,.2843 
.7761,.2846 
.7574,.2844 
.7403,.2840 
.7247,.2834 
.7103,.2825 
.6969,.2814 
.6846,.2801 
.6732,.2787 
.6625,.2770 
.6545,.2756 


m = —4 
.7966,.2843 
.7692,.2847 
.7442,.2847 
.7214,.2845 
.7003,.2840 
.6809,.2834 
.6628,.2826 
.6459,.2817 
.6300,.2807 
.6152,.2796 
.6039,.2787 


m = 2 
.7727,.4799 
.7890,.4788 
.8074,.4767 
.8283,.4735 
.8524,.4687 
.8802,.4618 
.9131,.4518 
.9526,.4367 
1.0020,.4123 
1.0667,.3665 

1.1283,.2730 


; = 4, n = 2 


m = 1 
.7727,.4799 
.7812,.4791 
.7913,.4775 
.8033,.4749 
.8176,.4710 
.8344,.4654 
.8543,.4573 
.8780,.4456 
.9065,.4276 
.9402,.3972 
.9655,.3535 


m = 0 
.7727,.4799 
.7734,.4795 
.7756,.4784 
.7794,.4763 
.7847,.4733 
.7918,.4690 
.8008,.4630 
.8119,.4546 
.8253,.4426 
.8404,.4247 
.8520,.4037 


m = —1 
.7727,.4799 
.7658,.4799 
.7604,.4792 
.7563,.4778 
.7536,.4756 
.7522,.4725 
.7522,.4683 
.7534,.4627 
.7558,.4554 
.7591,.4455 
.7619,.4352 


m = —2 
.7727,.4799 
.7583,.4803 
.7455,.4800 
.7342,.4792 
.7243,.4777 
.7155,.4757 
.7078,.4730 
.7012,.4696 
.6954,.4654 
.6904,.4602 
.6868,.4552 


m = —3 
.7727,.4799 
.7509,.4807 
.7311,.4808 
.7130,.4805 
.6965,.4797 
.6813,.4785 
.6672,.4769 
.6542,.4750 
.6422,.4727 
.6310,.4701 
.6225,.4678 


m = —4 
.7727,.4799 
.7436,.4811 
.7170,.4816 
.6927,.4817 
.6702,.4814 
.6493,.4808 
.6299,.4800 
.6117,.4789 
.5946,.4777 
.5786,.4764 
.5664,.4752 


The first few quasi-normal frequencies for extremal black holes asymptote to m/2 for m > 0, but imposing this 
behavior (for instance by assuming fitting functions of the form T ~ m/ 2 -|-/ 2 (l —usually decreases the accuracy 
of the fits. Indeed our ht for Z = 2, m = 2 performs better than Echeverria’s, and the one for Z = 2, m = 0 performs 
better than the fit provided in j^. This is shown in Fig. I21L where we plot the percentage error for each of the fits, 
in the range 0 < j < 0.99. 

For the I = \m\ modes we can also make contact with an approach initiated by Press and Goebel [b^, and 
investigated in depth by Mashhoon [^, in which QNMs are regarded as waves orbiting around the unstable photon 
orbit, and slowly leaking out. Mashhoon’s analytical formulas were recently compared to numerical results |64l| and 
found to be in good agreement for I = \m\. This approach basically predicts that J-i±io = \m\Muj n, w here loq is 
the co-rotating (counter-rotating) orbital frequency at the photon orbit for positive (negative) m (see for further 
details). The counter-rotating frequency is well approximated by Mujo ^ -1/7-1- 5{j — 1)/147. This implies that the 
fundamentaU = 2,m = —2 QNM frequency should be lF 2 _ 2 o = 2/7-1-10(1—j)/147 ~ 0.286-1-0.068(1—/). To compare 
with this formula, we can do a slightly different ht to the data for m = —2. If we fit the data to T 2 - 2 Q = b + c{l — j) we 
get b = 0.287, c = 0.0805 with a maximum error of 1.6%. This is in reasonable agreement with Mashhoon’s prediction. 
For / = 3, m = —3 the same fit yields T 3-30 = 0.456 -I- 0.134(1 — j) with 1.6% maximum error, whereas Mashhoon’s 
prediction is .^ 3-30 = 0.428 -I- 0.102(1 — j). Thus for negative m, I = |m| we can approximate the frequencies by 
J-i-io ~ \m\ [1/7 -I- 5(1 — j)/147]. For positive m we don’t get such good agreement with a linear ht, because Mui is 
a rapidly varying function of j. Nevertheless, for moderate to large j and I = m we get, using Mashhoon’s formula. 
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TABLE V: Angular separation constants for the first three overtones with I — 2. 
Tables we omit leading zeros. 


j m = 2 
0.00 4.0000,.0000 
0.10 3.8957,.0242 
0.20 3.7810,.0492 
0.30 3.6531,.0751 
0.40 3.5087,.1019 
0.50 3.3423,.1292 
0.60 3.1454,.1567 
0.70 2.9032,.1832 
0.80 2.5853,.2053 
0.90 2.1098,.2111 
0.98 1.3336,.1500 


I — 2, n — 0 

m = 1 m — 0 m = —1 

4.0000,.0000 4.0000,.0000 4.0000,.0000 

3.9485.. 0122 3.9993,.0003 4.0483,-.0115 

3.8932.. 0252 3.9972,.0014 4.0939,-.0222 

3.8332.. 0389 3.9937,.0031 4.1371,-.0322 

3.7676.. 0532 3.9886,.0056 4.1784,-.0415 

3.6947.. 0682 3.9817,.0088 4.2178,-.0500 

3.6125.. 0835 3.9730,.0126 4.2558,-.0577 

3.5172.. 0986 3.9619,.0172 4.2924,-.0645 

3.4023.. 1122 3.9480,.0223 4.3279,-.0705 

3.2534.. 1195 3.9304,.0276 4.3623,-.0755 
3.0797,.1022 3.9127,.0315 4.3892,-.0787 


j m = 2 
0.00 4.0000,.0000 
0.10 3.9031,.0741 
0.20 3.7958,.1504 
0.30 3.6756,.2287 
0.40 3.5386,.3091 
0.50 3.3791,.3908 
0.60 3.1882,.4726 
0.70 2.9498,.5513 
0.80 2.6313,.6171 
0.90 2.1455,.6342 
0.98 1.3457,.4499 


_ 1^2, n^l _ 

m = l m — 0 m = —1 

4.0000,.0000 4.0000,.0000 4.0000,.0000 

3.9524.. 0375 3.9998,.0010 4.0451,-.0355 

3.9017.. 0771 3.9990,.0040 4.0883,-.0690 

3.8469.. 1185 3.9977,.0090 4.1300,-.1004 

3.7869.. 1617 3.9956,.0161 4.1707,-.1296 

3.7197.. 2065 3.9925,.0252 4.2105,-.1566 

3.6429.. 2520 3.9880,.0364 4.2497,-.1813 

3.5523.. 2968 3.9815,.0496 4.2884,-.2034 

3.4405.. 3365 3.9722,.0645 4.3266,-.2227 

3.2918.. 3564 3.9587,.0800 4.3642,-.2389 

3.1248.. 2952 3.9442,.0906 4.3937,-.2494 


To save space, in this and the following 


m = —2 
4.0000,.0000 
4.0956,-.0233 
4.1839,-.0457 
4.2659,-.0674 
4.3427,-.0882 
4.4147,-.1083 
4.4827,-.1277 
4.5470,-.1463 
4.6081,-.1643 
4.6662,-.1816 
4.7108,-.1950 


m = —2 
4.0000,.0000 
4.0885,-.0720 
4.1702,-.1420 
4.2461,-.2100 
4.3172,-.2761 
4.3841,-.3405 
4.4476,-.4031 
4.5080,-.4641 
4.5658,-.5236 
4.6213,-.5815 
4.6643,-.6268 


Z = 2, n = 2 


J 


m = 2 


m = 1 


m — 0 


m = —1 


m = —2 


0.00 

0.10 

0.20 

0.30 

0.40 

0.50 

0.60 

0.70 

0.80 

0.90 

0.98 


. 0000,.0000 
.9158,.1285 
.8222,.2589 
.7162,.3912 
.5935,.5253 
.4480,.6602 
.2697,.7944 
.0405,.9229 
7230,1.0312 
2175,1.0596 
.3696,.7498 


4.0000, 

3.9594, 

3.9171, 

3.8721, 

3.8226, 

3.7666, 

3.7006, 

3.6195, 

3.5141, 

3.3653, 

3.1982, 


.0000 

0651 

1329 

.2031 

.2756 

.3498 

.4246 

4972 

5602 

.5869 

.4474 


4.0000, 

4.0007, 

4.0028, 

4.0061, 

4.0104, 

4.0150, 

4.0193, 

4.0222, 

4.0219, 

4.0163, 

4.0087, 


0000 

0015 

0061 

0137 

0246 

0388 

0563 

0771 

1005 

1237 

1378 


4.0000 

4.0399,- 

4.0798,- 

4.1202,- 

4.1616,- 

4.2041,- 

4.2479,- 

4.2928,- 

4.3384,- 

4.3835,- 

4.4184,- 


.0000 

.0624 

.1219 

.1785 

.2319 

.2817 

.3276 

.3691 

.4056 

.4366 

.4576 


4.0000,.0000 

4.0770,-.1266 

4.1486,-.2516 

4.2160,-.3750 

4.2805,-.4970 

4.3428,-.6178 

4.4039,-.7375 

4.4643,-.8560 

4.5245,-.9736 

4.5851,-1.0903 

4.6340,-1.1830 


Tm--\l\ (l/2-V3(l-j)/8). 


This is in very good agreement with the numerical data. 
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TABLE VI: Angular separation constants for the first three overtones with 1 = 3. 

I = 3, n = 0 

j m = 3 m = 2 m = 1 m = 0 m = —1 m = —2 m = —3 

0.00 1.0000,.0000 1.0000,.0000 1.0000,.0000 1.0000,.0000 1.0000,.0000 1.0000,.0000 1.0000,.0000 

0.10 9.8739,.0191 9.9167,.0128 9.9583,.0066 9.9988,.0004 1.0383,-.0058 1.0769,-.0119 1.1144,-.0180 

0.20 9.7338,.0392 9.8257,.0265 9.9128,.0139 9.9953,.0015 1.0736,-.0108 1.1481,-.0230 1.2190,-.0351 

0.30 9.5764,.0605 9.7257,.0410 9.8629,.0220 9.9893,.0033 1.1062,-.0150 1.2145,-.0332 1.3153,-.0512 

0.40 9.3970,.0828 9.6146,.0563 9.8080,.0307 9.9808,.0059 1.1361,-.0185 1.2766,-.0426 1.4043,-.0665 

0.50 9.1890,.1060 9.4895,.0722 9.7470,.0401 9.9695,.0092 1.1637,-.0211 1.3348,-.0511 1.4870,-.0809 

0.60 8.9414,.1299 9.3461,.0884 9.6785,.0500 9.9552,.0132 1.1890,-.0229 1.3896,-.0587 1.5642,-.0946 

0.70 8.6358,.1534 9.1773,.1043 9.6006,.0601 9.9375,.0179 1.2121,-.0238 1.4412,-.0656 1.6364,-.1075 

0.80 8.2345,.1739 8.9702,.1181 9.5096,.0697 9.9157,.0231 1.2330,-.0239 1.4900,-.0716 1.7043,-.1198 

0.90 7.6364,.1815 8.6940,.1235 9.3992,.0767 9.8889,.0285 1.2518,-.0231 1.5362,-.0768 1.7682,-.1314 

0.98 6.6687,.1319 8.3443,.0951 9.2872,.0749 9.8630,.0326 1.2653,-.0218 1.5713,-.0803 1.8167,-.1403 

I = 3, n = 1 

j m = 3 m = 2 m = 1 m = 0 m = —1 m = —2 m = —3 

0.00 1.0000,.0000 1.0000,.0000 1.0000,.0000 1.0000,.0000 1.0000,.0000 1.0000,.0000 1.0000,.0000 

0.10 9.8775,.0578 9.9192,.0388 9.9597,.0199 9.9991,.0011 1.0375,-.0176 1.0749,-.0363 1.1113,-.0548 

0.20 9.7413,.1186 9.8312,.0800 9.9161,.0420 9.9965,.0044 1.0726,-.0329 1.1448,-.0700 1.2132,-.1069 

0.30 9.5881,.1825 9.7344,.1237 9.8686,.0663 9.9920,.0098 1.1056,-.0460 1.2104,-.1013 1.3072,-.1563 

0.40 9.4131,.2495 9.6268,.1697 9.8165,.0927 9.9854,.0175 1.1366,-.0566 1.2723,-.1301 1.3944,-.2033 

0.50 9.2094,.3192 9.5053,.2174 9.7585,.1209 9.9766,.0273 1.1658,-.0649 1.3310,-.1565 1.4758,-.2480 

0.60 8.9659,.3905 9.3652,.2661 9.6931,.1505 9.9651,.0392 1.1934,-.0707 1.3870,-.1805 1.5522,-.2905 

0.70 8.6633,.4607 9.1992,.3136 9.6180,.1807 9.9504,.0531 1.2193,-.0740 1.4405,-.2020 1.6240,-.3310 

0.80 8.2629,.5217 8.9935,.3544 9.5294,.2090 9.9316,.0685 1.2435,-.0747 1.4918,-.2212 1.6919,-.3696 

0.90 7.6604,.5443 8.7162,.3698 9.4211,.2293 9.9077,.0846 1.2658,-.0726 1.5410,-.2378 1.7563,-.4064 

0.98 6.6779,.3956 8.3645,.2808 9.3156,.2228 9.8841,.0963 1.2820,-.0690 1.5790,-.2494 1.8055,-.4346 

I = 3, n = 2 

j m = 3 m = 2 m = 1 m = 0 m = —1 m = —2 m = —3 

0.00 1.0000,.0000 1.0000,.0000 1.0000,.0000 1.0000,.0000 1.0000,.0000 1.0000,.0000 1.0000,.0000 

0.10 9.8842,.0980 9.9238,.0658 9.9623,.0337 9.9997,.0018 1.0361,-.0301 1.0713,-.0619 1.1056,-.0937 

0.20 9.7555,.2006 9.8415,.1355 9.9225,.0710 9.9989,.0071 1.0710,-.0566 1.1389,-.1200 1.2028,-.1833 

0.30 9.6105,.3078 9.7512,.2089 9.8798,.1118 9.9974,.0159 1.1051,-.0794 1.2035,-.1743 1.2933,-.2692 

0.40 9.4442,.4196 9.6505,.2857 9.8331,.1559 9.9950,.0283 1.1385,-.0984 1.2657,-.2248 1.3780,-.3515 

0.50 9.2492,.5353 9.5360,.3653 9.7811,.2030 9.9911,.0443 1.1715,-.1134 1.3262,-.2716 1.4580,-.4305 

0.60 9.0138,.6534 9.4027,.4460 9.7218,.2522 9.9853,.0637 1.2040,-.1245 1.3853,-.3146 1.5341,-.5064 

0.70 8.7175,.7691 9.2422,.5244 9.6524,.3020 9.9767,.0864 1.2359,-.1314 1.4432,-.3538 1.6068,-.5794 

0.80 8.3192,.8696 9.0394,.5910 9.5687,.3485 9.9641,.1116 1.2667,-.1339 1.5002,-.3892 1.6766,-.6498 

0.90 7.7083,.9064 8.7596,.6140 9.4645,.3796 9.9459,.1375 1.2960,-.1319 1.5562,-.4207 1.7439,-.7176 

0.98 6.6964,.6591 8.3986,.4546 9.3745,.3656 9.9274,.1558 1.3177,-.1274 1.6004,-.4431 1.7963,-.7701 




j j 


FIG. 21: Accnracy of the fits of frequency (black) and quality factor (red) for the fundamental mode with I = m = 2 (left) and 
I = 2, m = 0 (right). Solid lines refer to our fits, dashed lines to the fits proposed in Refs. |3 and 1^. 
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TABLE VII: Angular separation constants for the first three overtones with I = 4. 


t = 4, n = 0 


j m = 4 m = 3 m = 2 

0.00 18.0000,.0000 18.0000,.0000 18.0000,.0000 
0.10 17.8633,.0156 17.8978,.0118 17.9317,.0081 
0.20 17.7101,.0323 17.7842,.0248 17.8555,.0172 
0.30 17.5365,.0502 17.6569,.0388 17.7703,.0275 
0.40 17.3369,.0694 17.5128,.0539 17.6743,.0387 
0.50 17.1033,.0897 17.3475,.0698 17.5655,.0509 
0.60 16.8229,.1110 17.1544,.0865 17.4405,.0637 
0.70 16.4735,.1324 16.9225,.1029 17.2943,.0768 
0.80 16.0106,.1519 16.6313,.1174 17.1183,.0888 
0.90 15.3146,.1609 16.2310,.1230 16.8945,.0959 
0.98 14.1819,.1190 15.6920,.9097 16.6463,.0850 


m = l m = 0 m = —1 

18.0000,.0000 18.0000,.0000 18.0000,.0000 

17.9650.. 0043 17.9977,.0006 18.0298,-.0032 

17.9243.. 0097 17.9906,.0022 18.0547,-.0053 

17.8774.. 0162 17.9787,.0049 18.0749,-.0064 

17.8235.. 0237 17.9618,.0088 18.0905,-.0063 

17.7616.. 0323 17.9394,.0137 18.1017,-.0051 

17.6903.. 0417 17.9111,.0197 18.1083,-.0028 

17.6075.. 0518 17.8762,.0266 18.1104,.0006 

17.5099.. 0619 17.8335,.0343 18.1076,.0050 

17.3918.. 0707 17.7814,.0425 18.0997,.0105 

17.2747.. 0734 17.7313,.0487 18.0894,.0155 


m = —2 
18.0000,.0000 
18.0614,-.0070 
18.1167,-.0129 
18.1663,-.0177 
18.2109,-.0215 
18.2508,-.0242 
18.2862,-.0259 
18.3175,-.0265 
18.3448,-.0261 
18.3682,-.0247 
18.3843,-.0228 


m = —3 
18.0000,.0000 
18.0925,-.0108 
18.1766,-.0205 
18.2535,-.0292 
18.3239,-.0370 
18.3885,-.0438 
18.4480,-.0497 
18.5029,-.0547 
18.5536,-.0588 
18.6003,-.0622 
18.6352,-.0642 


m = —4 
18.0000,.0000 
18.1231,-.0146 
18.2347,-.0281 
18.3366,-.0408 
18.4302,-.0527 
18.5164,-.0637 
18.5963,-.0741 
18.6706,-.0838 
18.7399,-.0929 
18.8046,-.1013 
18.8535,-.1078 


i = 4, n = 1 


j 

m = 4 

m = 3 

m = 

2 

m = 

1 

m — 0 

m = - 

-1 

0.00 

18.0000,.0000 

18.0000,.0000 

18.0000, 

.0000 

18.0000, 

.0000 

18.0000,.0000 

18.0000,. 

0000 

0.10 

17.8655,.0471 

17.8995,.0357 

17.9330, 

.0244 

17.9658, 

.0130 

17.9980,.0016 

18.0296,- 

.0097 

0.20 

17.7148,.0974 

17.7881,.0746 

17.8586, 

.0519 

17.9265, 

.0293 

17.9919,.0065 

18.0550,- 

.0162 

0.30 

17.5441,.1512 

17.6633,.1168 

17.7756, 

.0827 

17.8816, 

.0487 

17.9817,.0147 

18.0764,- 

.0194 

0.40 

17.3476,.2086 

17.5220,.1620 

17.6823, 

.1164 

17.8301, 

.0713 

17.9669,.0261 

18.0939,- 

.0194 

0.50 

17.1173,.2695 

17.3596,.2100 

17.5762, 

.1529 

17.7710, 

.0969 

17.9473,.0408 

18.1075,- 

.0160 

0.60 

16.8400,.3331 

17.1693,.2597 

17.4539, 

.1914 

17.7027, 

.1251 

17.9221,.0586 

18.1172,- 

.0094 

0.70 

16.4932,.3974 

16.9396,.3091 

17.3102, 

.2305 

17.6228, 

.1552 

17.8906,.0792 

18.1229,. 

0005 

0.80 

16.0315,.4557 

16.6493,.3522 

17.1358, 

.2663 

17.5279, 

.1856 

17.8514,.1023 

18.1241,. 

0135 

0.90 

15.3329,.4824 

16.2473,.3686 

16.9127, 

.2872 

17.4121, 

.2115 

17.8028,.1265 

18.1204,. 

0296 

0.98 

14.1893,.3570 

15.7029,.2713 

16.6706, 

.2524 

17.2991, 

.2191 

17.7555,.1447 

18.1135,. 

0443 


m = —2 
18.0000,.0000 
18.0607,-.0211 
18.1159,-.0391 
18.1662,-.0539 
18.2120,-.0655 
18.2538,-.0740 
18.2919,-.0794 
18.3264,-.0817 
18.3575,-.0809 
18.3853,-.0771 
18.4051,-.0719 


m = —3 
18.0000,.0000 
18.0913,-.0326 
18.1747,-.0621 
18.2514,-.0886 
18.3222,-.1122 
18.3879,-.1332 
18.4490,-.1514 
18.5061,-.1671 
18.5595,-.1802 
18.6095,-.1910 
18.6474,-.1979 


m = —4 
18.0000,.0000 
18.1212,-.0440 
18.2314,-.0852 
18.3323,-.1237 
18.4251,-.1598 
18.5110,-.1936 
18.5908,-.2254 
18.6654,-.2552 
18.7352,-.2833 
18.8009,-.3097 
18.8507,-.3296 


i = 4, n = 2 


j 

m = 4 

m = 3 

m = 

2 

m = 

1 

m — 0 

m = - 

-1 

0.00 

18.0000,.0000 

18.0000,.0000 

18.0000, 

.0000 

18.0000, 

.0000 

18.0000,.0000 

18.0000,. 

0000 

0.10 

17.8697,.0793 

17.9029,.0601 

17.9355, 

.0410 

17.9674, 

.0219 

17.9987,.0027 

18.0293,- 

.0165 

0.20 

17.7240,.1637 

17.7958,.1255 

17.8647, 

.0873 

17.9309, 

.0490 

17.9946,.0107 

18.0558,- 

.0277 

0.30 

17.5589,.2537 

17.6760,.1960 

17.7862, 

.1387 

17.8899, 

.0815 

17.9876,.0241 

18.0796,- 

.0336 

0.40 

17.3686,.3494 

17.5403,.2715 

17.6981, 

.1950 

17.8434, 

.1191 

17.9773,.0429 

18.1009,- 

.0341 

0.50 

17.1448,.4505 

17.3836,.3513 

17.5976, 

.2559 

17.7899, 

.1616 

17.9632,.0669 

18.1196,- 

.0293 

0.60 

16.8739,.5561 

17.1987,.4341 

17.4808, 

.3200 

17.7275, 

.2086 

17.9444,.0962 

18.1356,- 

.0191 

0.70 

16.5324,.6625 

16.9734,.5159 

17.3418, 

.3849 

17.6535, 

.2586 

17.9198,.1302 

18.1485,- 

.0034 

0.80 

16.0730,.7590 

16.6851,.5871 

17.1708, 

.4440 

17.5639, 

.3087 

17.8876,.1681 

18.1578,. 

0174 

0.90 

15.3693,.8032 

16.2796,.6133 

16.9490, 

.4770 

17.4528, 

.3507 

17.8458,.2076 

18.1626,. 

0431 

0.98 

14.2040,.5948 

15.7228,.4476 

16.7200, 

.4106 

17.3486, 

.3620 

17.8045,.2369 

18.1625,. 

0666 


m = —2 
18.0000,.0000 
18.0594,-.0358 
18.1146,-.0664 
18.1662,-.0919 
18.2148,-.1123 
18.2607,-.1276 
18.3042,-.1379 
18.3454,-.1432 
18.3844,-.1435 
18.4210,-.1389 
18.4484,-.1317 


m = —3 
18.0000,.0000 
18.0889,-.0551 
18.1711,-.1053 
18.2477,-.1507 
18.3197,-.1916 
18.3877,-.2281 
18.4524,-.2604 
18.5142,-.2886 
18.5735,-.3129 
18.6306,-.3333 
18.6748,-.3469 


m = —4 
18.0000,.0000 
18.1178,-.0744 
18.2254,-.1444 
18.3244,-.2102 
18.4161,-.2722 
18.5017,-.3308 
18.5820,-.3861 
18.6577,-.4383 
18.7293,-.4879 
18.7975,-.5348 
18.8498,-.5707 


TABLE VIII: Eitting coefficients for the dimensionless frequency, Timn = /i + /2(1 — , and the quality factor Qimn = 

qi +?2(1 —i)® of a Kerr black hole. We give coefficients for I = 2, all values of m, and the three lowest overtones (n = 0, 1, 2). 
Eor each fit we also give the maximum percentage error in the range j € [0, 0.99]. 


m 

n 

/l /2 

/3 

% 


?2 

<13 

% 

2 

0 

1.5251 -1.1568 

0.1292 

1.85 

0.7000 

1.4187 

-0.4990 

0.88 


1 

1.3673 -1.0260 

0.1628 

1.56 

0.1000 

0.5436 

-0.4731 

1.69 


2 

1.3223 -1.0257 0.1860 

1.91 

-0.1000 

0.4206 

-0.4256 

2.52 

1 

0 

0.6000 -0.2339 

0.4175 

2.03 

-0.3000 

2.3561 

-0.2277 

3.65 


1 

0.5800 -0.2416 

0.4708 

2.40 

-0.3300 

0.9501 

-0.2072 

3.18 


2 

0.5660 -0.2740 

0.4960 

4.04 

-0.1000 

0.4173 

-0.2774 

2.46 

0 

0 

0.4437 -0.0739 

0.3350 

1.04 

4.0000 

-1.9550 

0.1420 

2.63 


1 

0.4185 -0.0768 

0.4355 

1.50 

1.2500 

-0.6359 

0.1614 

4.01 


2 

0.3734 -0.0794 0.6306 

2.72 

0.5600 

-0.2589 

0.3034 

4.33 

-1 

0 

0.3441 0.0293 

2.0010 

0.07 

2.0000 

0.1078 

5.0069 

2.82 


1 

0.3165 0.0301 

2.3415 

0.05 

0.6100 

0.0276 

13.1683 

0.67 


2 

0.2696 0.0315 

2.7755 

0.43 

0.2900 

0.0276 

6.4715 

2.40 

-2 

0 

0.2938 0.0782 

1.3546 

0.63 

1.6700 

0.4192 

1.4700 

0.71 


1 

0.2528 0.0921 

1.3344 

0.87 

0.4550 

0.1729 

1.3617 

0.79 


2 

0.1873 0.1117 

1.3322 

1.34 

0.1850 

0.1266 

1.3661 

1.16 
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TABLE IX: Same as Table IVmi but for I = 3. 


m 

n 

/i 

/2 

/3 

% 

<11 

<72 

<13 

% 

3 

0 

1.8956 -1.3043 

0.1818 

1.36 

0.9000 

2.3430 

-0.4810 

0.42 


1 

1.8566 -1.2818 

0.1934 

1.35 

0.2274 

0.8173 

-0.4731 

0.88 


2 

1.8004 -1.2558 

0.2133 

1.28 

0.0400 

0.5445 

-0.4539 

1.52 

2 

0 

1.1481 -0.5552 

0.3002 

1.09 

0.8313 

2.3773 

-0.3655 

1.28 


1 

1.1226 -0.5471 

0.3264 

1.23 

0.2300 

0.8025 

-0.3684 

0.51 


2 

1.0989 -0.5550 

0.3569 

1.41 

0.1000 

0.4804 

-0.3784 

0.81 

1 

0 

0.8345 -0.2405 

0.4095 

1.12 

23.8450 

-20.7240 

0.03837 

3.47 


1 

0.8105 -0.2342 

0.4660 

1.55 

8.8530 

-7.8506 

0.03418 

3.64 


2 

0.7684 -0.2252 

0.5805 

2.67 

2.1800 

-1.6273 

0.1136 

4.04 

0 

0 

0.6873 -0.09282 

0.3479 

0.83 

6.7841 

-3.6112 

0.09480 

3.99 


1 

0.6687 -0.09155 

0.4021 

0.95 

2.0075 

-0.9930 

0.12297 

4.18 


2 

0.6343 -0.08915 

0.5117 

1.28 

0.9000 

-0.3409 

0.2679 

2.89 

-1 

0 

0.5751 0.02508 

3.1360 

0.42 

3.0464 

0.1162 

-0.2812 

2.65 


1 

0.5584 0.02514 

3.4154 

0.42 

1.2000 

-0.1928 

0.1037 

2.75 


2 

0.5271 0.02561 

3.8011 

0.29 

1.0000 

-0.4424 

0.02467 

3.15 

-2 

0 

0.5158 0.08195 

1.4084 

0.35 

2.9000 

0.3356 

2.3050 

0.72 


1 

0.4951 0.08577 

1.4269 

0.41 

0.9000 

0.1295 

1.6142 

0.80 


2 

0.4567 0.09300 

1.4469 

0.53 

0.4900 

0.0848 

1.9737 

0.52 

-3 

0 

0.4673 

0.1296 

1.3255 

0.61 

2.5500 

0.6576 

1.3378 

0.79 


1 

0.4413 

0.1387 

1.3178 

0.68 

0.7900 

0.2381 

1.3706 

0.73 


2 

0.3933 

0.1555 

1.3037 

0.82 

0.4070 

0.1637 

1.3819 

0.88 


TABLE X: Same as Table IVmi but for I — 4. 


m 

n 

/l /2 

/3 

% 

<11 

<12 

<13 

% 

4 

0 

2.3000 -1.5056 

0.2244 

1.83 

1.1929 

3.1191 

-0.4825 

0.37 


1 

2.3000 -1.5173 

0.2271 

1.75 

0.3000 

1.1034 

-0.4703 

0.77 


2 

2.3000 -1.5397 0.2321 

1.61 

0.1100 

0.6997 

-0.4607 

0.10 

3 

0 

1.6869 -0.8862 

0.2822 

1.05 

1.4812 

2.8096 

-0.4271 

0.14 


1 

1.6722 -0.8843 

0.2923 

1.10 

0.4451 

0.9569 

-0.4250 

0.37 


2 

1.6526 -0.8888 

0.3081 

1.15 

0.2200 

0.5904 

-0.4236 

0.66 

2 

0 

1.2702 -0.4685 

0.3835 

1.11 

-3.6000 

7.7749 

-0.1491 

0.97 


1 

1.2462 -0.4580 

0.4139 

1.39 

-1.5000 

2.8601 

-0.1392 

0.12 


2 

1.2025 -0.4401 

0.4769 

2.26 

-1.5000 

2.2784 

-0.1124 

0.31 

1 

0 

1.0507 -0.2478 

0.4348 

0.97 

14.0000 

-9.8240 

0.09047 

0.81 


1 

1.0337 -0.2439 

0.4695 

1.15 

4.2000 

-2.8399 

0.1081 

0.91 


2 

1.0019 -0.2374 0.5397 

1.53 

2.2000 

-1.4195 

0.1372 

0.53 

0 

0 

0.9175 -0.1144 0.3511 

0.75 

7.0000 

-2.7934 

0.1708 

0.26 


1 

0.9028 -0.1127 0.3843 

0.82 

2.2000 

-0.8308 

0.2023 

0.26 


2 

0.8751 -0.1096 

0.4516 

0.96 

1.2000 

-0.4159 

0.2687 

0.60 

-1 

0 

0.7908 0.02024 

5.4628 

0.96 

4.6000 

-0.4038 

0.4629 

2.52 


1 

0.7785 0.02005 

5.8547 

0.98 

1.6000 

-0.2323 

0.2306 

2.37 


2 

0.7549 0.01985 

6.5272 

0.96 

1.6000 

-0.8136 

0.03163 

2.32 

-2 

0 

0.7294 0.07842 

1.5646 

0.23 

4.0000 

0.2777 

2.0647 

2.11 


1 

0.7154 0.07979 

1.5852 

0.25 

1.3200 

0.08694 

4.3255 

0.75 


2 

0.6885 0.08259 

1.6136 

0.32 

0.7500 

0.05803 

3.7971 

0.66 

-3 

0 

0.6728 0.1338 

1.3413 

0.43 

3.700 

0.5829 

1.6681 

0.45 


1 

0.6562 0.1377 

1.3456 

0.46 

1.1800 

0.2111 

1.4129 

0.70 


2 

0.6244 0.1454 

1.3513 

0.52 

0.6600 

0.1385 

1.3742 

0.82 

-4 

0 

0.6256 0.1800 

1.3218 

0.62 

3.4000 

0.8696 

1.4074 

0.63 


1 

0.6061 0.1869 

1.3168 

0.67 

1.0800 

0.3095 

1.3279 

0.81 


2 

0.5686 0.2003 

1.3068 

0.74 

0.5980 

0.2015 

1.3765 

0.69 
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